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ABSTRACT 


This  report  examines  optimum  array  processing  for  bit  detection 
of  a  binary  communication  signal  in  the  presence  of  a  directional 
interfering  signal.  The  binary  communication  signal  is  assumed  to 
be  either  completely  known  or  to  have  random  parameters  (Rayleigh 
amplitude  and  uniform  phase).  The  Interfering  signal  is  a  wide- 
sense  stationary  Gaussian  process.  The  receiving  antenna  is  a  linear 
array  of  equally  spaced  isotropic  point  elements.  Statistically 
independent  white  noise  is  added  to  the  signal  at  each  element. 

An  analytical  expression  for  the  inverse  of  the  covariance 
matrix  is  obtained,  and  the  general  structure  for  the  optimum 
detector  is  derived.  The  optimum  processor  is  a  correlation 
receiver  whose  detailed  structure  is  dependent  in  general  on  the  power 
spectrum  of  the  interference,  as  well  as  the  other  antenna  and  signal 
parameters. 

The  detection  error  is  calculated  for  several  cases.  It  is 
found  that  the  detection  error  is  minimized  if  the  array  elements 
have  an  optimum  spacing,  determined  by  the  array  size  and  the  in¬ 
cidence  angles  of  the  desired  signal  and  the  interference.  When 
the  optimum  spacing  is  used,  the  structure  of  the  detection  system 
is  not  dependent  on  the  power  spectrum  of  the  interfering  signal. 
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CHAPTER  I 
INTRODUCTION 


Statistical  decision  theory  has  been  applied  for  some  time 
to  the  problem  of  determining  the  optimum  structure  of  a  receiver. 

The  basic  concept  of  optimal  detection  for  completely  known  signals 
or  known  signal  with  unknown  parameters  corrupted  by  additive 
Gaussian  noise  has  been  discussed  extensively  in  D-5].  During 
the  past  decade  this  theory  has  been  extended  to  analyze  a  com¬ 
munication  system  which  includes  an  antenna  array.  The  philosophy 
is  to  seek  an  optimum  performance  of  the  entire  system  based  on 
statistical  decision  theory,  rather  than  to  consider  the  performance 
of  the  array  or  the  receiver  separately. 

Basically,  a  decision  theory  receiver  processes  the  received 
signals  and  makes  one  of  several  possible  decisions,  based  on  an 
optimum  decision  rule.  Peterson  et  al  [63  showed  that  the  re¬ 
ceiver  which  forms  the  likelihood  ratio  from  the  received  signal 
and  compares  it  with  a  threshold  is  an  optimal  receiver  regardless 
of  which  criterion  is  used  (e.g.,  Bayes,  minimax,  Neyman-Pearson, 
etc.).  The  various  criteria  affect  the  receiver  only  in  the  value 
of  the  threshold  setting.  A  paper  by  Stocklin  C73  formally  showed 
that  optimal  space-time  signal  processing  also  called  for  likelihood 
ratio  detection.  Bryn  C81  applied  Rice's  Fourier  series  repre¬ 
sentation  to  a  vector  representation  of  the  signals  to  obtain  a 
likelihood  ratio  detection  good  for  wi de-sense  stationary  processes. 
The  problem  of  spacing  a  finite  number  of  point  detectors  in  a  plane 
to  detect  spatially  isotropic  signals  corrupted  by  spatially  iso¬ 
tropic  and  covariance-separable  noise  was  studied  extensively  by 
Gaerder  C9J.  Capon  CIO],  Young  Cll-13],  Gallop  C14]  have  also  treated 
various  problems  in  this  area. 

The  goal  of  this  research  is  to  find  and  examine  the  behavior 
of  an  optimum  detector  for  detecting  binary  signals  received  by  a 
linear  array  of  equally  spaced  isotropic  point  elements.  The  novelty 
of  the  problem  treated  here  is  that  the  signal  is  assumed  to  be 
corrupted  by  an  external  directional  interfering  signal.  The 
interfering  signal  is  modeled  as  a  wide-sense  stationary  process. 

In  addition,  the  processing  unit  behind  each  antenna  is  assumed  to 
contribute  white  noise,  statistically  independent  from  one  element 
to  the  next. 

The  research  discussed  here  is  motivated  by  the  adaptive 
antenna  array  research  currently  being  conducted  at  the  ElectroScience 
LaboratoryC24-26].  The  purpose  of  this  research  is  to  determine 
the  performance  of  the  optimum  detection  system  so  that  this 
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performance  may  be  taken  as  a  base  for  evaluating  the  performance 
of  the  adaptive  antenna  array. 

A  similar  problem  has  been  treated  by  Urkowitz  HI 5].  However, 
his  approach  requires  the  solution  of  a  matrix  integral  equation, 
and  his  results,  which  are  carried  out  only  in  general  terms,  do 
not  include  a  detailed  examination  of  specific  cases. 

Two  different  series  representations  of  a  vector  random 
process  are  discussed  in  Van  Trees  C3],  These  two  methods  both 
require  the  solution  to  a  matrix  integral  equation,  which  is  tedious 
and  difficult.  The  method  used  here  is  the  scalar  Karhunen-Loeve 
representation.  The  advantage  of  this  choice  is  that  only  one 
integral  equation  needs  to  be  solved.  The  material  will  be  pre¬ 
sented  in  the  following  manner. 

Chapter  II  is  devoted  to  formulating  a  workable  expression  for  the 
covariance  matrix  and  its  inverse.  First,  the  inverse  of  the 
covariance  matrix  associated  with  broadside  interference  is  obtained 
by  generalizing  the  results  for  a  two-element  array,  a  three- 
element  array,  etc.  Then  an  expression  for  the  inverse  valid  for 
an  arbitrary  arrival  angle  is  obtained  by  employing  a  trans¬ 
formation.  This  transformation  relates  the  original  observation 
space  to  a  modified  space  associated  with  delayed  versions  of  the 
signals. 

In  Chapter  III,  the  structure  of  the  optimum  processor  for 
detecting  completely  known  binary  signals  is  found  from  the  likelihood 
ratio  test.  Also,  the  equations  for  the  receiver  operating  character¬ 
istics  are  given.  To  make  the  results  more  specific,  some  numerical 
results  for  the  case  when  the  desired  communication  signals  are 
biphase  modulated  are  given  in  Chapter  IV.  The  discussion  con¬ 
centrates  on  the  directional  characteristics  of  the  optimum  detector 
and  its  error  rate  performance. 

The  optimum  system  for  detecting  the  signal  which  has  a  random 
amplitude  and  phase  is  Investigated  in  Chapter  V.  We  assume  that  the 
amplitude  Is  a  Rayleigh  random  variable  and  the  phase  is  uniformly 
distributed  between  and  ir.  Some  numerical  results  for  the  case 
when  the  desired  signals  are  on-off  keying  signals  are  discussed  in 
detail  in  the  last  section  of  this  Chapter. 

In  Chapter  VI,  we  briefly  discuss  the  general  optimum  system  for 
detecting  binary  signals  with  phase  and  amplitude  distributed  arbi¬ 
trarily.  Also,  we  prove  that  the  processing  outputs  are  free  from 
Interference  when  an  optimum-spaced  array  Is  used  In  the  detection 
system. 
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CHAPTER  II 

THE  COVARIANCE  MATRIX  AND  ITS  INVERSE 


A.  Introduction 


The  mathematical  model  of  a  general  array  detection  system 
is  shown  in  Fig.  2-1.  The  antenna  array  which  is  linear,  and 
consists  of  m  equally  spaced  isotropic  point  elements,  is  designed 
to  receive  the  signal  x'(t)  from  a  remote  station.  The  signal 
x'(t)  may  be  either  one  of  two  waveforms, 

Ja'(t)  under  hypothesis  Ha 
X'(t)  =•<  a 

lb'(t)  under  hypothesis  Hb 

and  the  task  of  the  receiver  is  to  decide  which  waveform  was  sent. 
The  signal  transmitted  by  the  source  is  corrupted  by  an  additive 
directional  interference  n'(t).  The  output  at  each  antenna  terminal 
will  be  the  sum  of  the  desired  communication  signal  and  the  inter¬ 
fering  random  noise  component.  Those  two  components  are  temporal 
and  spatial  functions. 

Suppose  the  incident  angles  for  the  desired  and  interfering 
signals  measured  from  the  array  axis  are  e  and  <j >,  respectively. 

Then  the  signals  received  by  two  adjacent  antennas  will  be  out  of 
phase  due  to  differences  in  the  arrival  times.  The  time  delays 
t  and  v,  given  by 

(2-1 )  x  =  —  •  cos  0 


(2-2)  v  =  £  cos  <|»  , 

are  the  arrival  time  differences  between  adjacent  elements  for  the 
desired  signal  and  the  interference,  respectively,  d  is  the  spacing 
between  antennas,  and  c  is  the  speed  of  light. 

Often,  the  received  signal  is  amplified  to  a  certain  level 
before  being  processed.  We  assume  that  the  gain  and  phase  shift  for 
each  of  the  m  ideal,  linear  amplifiers  are  identical.  For  simplicity, 
but  without  losing  generality,  the  phase  shift  of  each  amplifier  is 
assumed  to  be  zero.  Let  x-j(t)  and  n-j(t)  be  the  outputs  of  the  ith 
amplifier  when  the  inputs  are  xj(t)  and  nj(t),  respectively.  Since 


the  amplifiers  are  ideal  and  do  not  introduce  phase  shift,  we  have 


x.(t)  =  x[t-(i-l )tJ 
(2-3)  1 

n.j  (t)  =  n[t-(i-l)v3  ,  1  1  i  lm* 

The  interfering  signal  n(t)  is  assumed  to  be  a  wide-sense 
stationary  Gaussian  process  of  zero-mean,*  i.e.. 


(2-4)  E(n(t) }  =  E(n.(t) }  =  0 


(2-5)  E(n(t)  •  n(z)}  i  R(t-z)  =  R(z-t) 

where  R(t-z)  is  the  autocorrelation  function  of  the  signal  n(t). 

The  functions  Wi(t),  1  <_  i  <  m,  in  Fig.  2-1  represent  internal 
noises,  generated  by  the  amplifiers  and  other  electronic  components 
(e.g.,  mixers)  in  each  processing  unit  (not  shown  in  the  figure). 

We  further  assume  that  the  interference  and  the  internal  noises 
are  statistically  independent,  and  wi(t)  is  a  Gaussian  white  noise 
process  with  zero-mean  and  spectral  height  hw  for  all  i.  Since 
the  noise  sources  in  each  channel  are  statistically  independent, 
we  have 


(2-6)  E{w. (t)  •  w j ( z ) }  =  •  hw  •  6(t-z) 

where  Sjj  is  the  Kronecker  delta  and  6(»)  is  the  Dirac  delta 
function. 

The  goal  of  this  research  is  to  find  an  optimum  processor  for 
detecting  the  signal  x(t)  in  the  presence  of  interference.  We 
treat  here  the  binary  detection  problem.  That  is,  the  signal  to 
be  processed  in  the  ith  channel,  rj(t),  is  of  one  of  two  forms, 
as  fol 1 ows : 

(2-7)  hypothesis  H_:  r.(t)  =  a.(t)  +  n.(t)  +  w.(t) 

a  1  1  1  I 

(2-8)  hypothesis  Hb:  r^(t)  =  b^ (t)  +  n..(t)  +  w^ (t) 


*E{ - }  is  the  expectation  operator. 
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B.  The  Covariance  Matrix 


There  are  many  ways  of  characterizing  waveforms  and  random 
processes.  „The  method  adopted  here  uses  an  orthogonal  series 
Karhunen-Loeve  representation  C17].  That  is,  we  express  the  wave¬ 
form  as  a  series  expansion  on  the  orthonormal  basis  function 
f 4>  (t) }  given  by  the  eigenfunctions  of  the  integral  equation  p8] 


(2-9) 


0£(u)  •  du. 


<  £  <  oo 


where  the  bounds  specify  the  observation  interval,  from  t  =  -T  to 
t=T,  and  the  kernel  of  the  integral  equation  is  given  by. 


(2-10)  c(t,u)  =  E{[n.(t)  +  w1(t)Ini(u)  +  w.(u)]} 

=  R(t-u)  +  hw  •  6(t-u),  1  <^i  <  m  . 


We  have  made  use  of  Eqs.  (2-5)  and  (2-6),  and  the  assumption  that 
n-j(t)  and  w,(t)  are  uncorrelated,  to  get  Eq.  (2-10).  The  integral 
equation  may  then  be  written 

(2-11)  (Vhw)  *  ♦*(*)  =  J  R(t-u)  **£(u)  .  du 

For  a  Gaussian  random  process,  the  coefficients  in  the  expansion  are 
statistically  independent  Gaussian  random  variables  [193.  It  is  in 
this  case  that  the  expansion  becomes  attractive. 

We  approximate  the  signal  r.(t)  for  the  ith  channel  by  an  N-term 
expansion  1 


(2-12)  rf(t)  =  ru  ♦„(!) 

N 

=  Ji (>u  +  "1i  +  w,jt) '  1  -1  -* 

under  hypothesis  H_,  or 

O 

N 

(2-13)  r^(t)  =  l  (bi£  +  n<£  +  w^£)  •  «j»£(t),  1  <  i  <  m 


".V  -l- 


p 

b 


under  hypothesis  Hb,  where 


(2-14) 


ri£ 

*r*.(t)‘ 

ai£ 

fT 

a.(t) 

bU 

M-T 

^(t) 

ni£ 

n.(t) 

."I*. 

wi(t)  _ 

^£(t)  •  dt 


We  define  the  observation  vector  R  for  the  detection  system  as 
Ri 


(2-15)  R  4 


where  Ri,  associated  with  the  ith  channel,  is  itself  a  column 
vector,  i.e.. 


ril 


i2 


(2-16)  Ri  a 


11 


LriNj 


* - K  1* - 2 - Jt.. 


The  dimension  for  the  Gaussian  random  vector  R  is  mN.  We  define* 


(2-17) 

M  =  E(R|Hk>  =  E 


Ri  i»k 
R2  lHk 


Rih 


>  A 


H, 


M„ 


Mi 


under  hypothesis  H. 

a 


[B  under  hypothesis  Hb 


and  the  covariance  matrix 


(2-18)  K  *  E  {(R-M)  (R-M)1}  a  Q"1 

If  K  is  nonsingular,**  the  probability  density  function  of 
an  m-channel,  N-term  expansion  is  C20] 

1  -  7  (R-M)TQ(R-M) 

(2-19)  p(R|H^^fT  'e 

where  |K|  is  the  determinant  of  K 


*  We  use  Hk  to  represent  either  hypothesis. 

**If  K  Is  singular,  we  have  a  singular  detection  problem, 
singular  If  the  spectral  height  h.  Is  zero.  We  will  not 
this  case  here.  " 


R  for 


K  becomes 
examine 


The  following  several  pages  will  be  devoted  to  finding  a  workable 
expression  for  the  covariance  matrix.  Symbolically,  we  may  write 


K  = 

(2-20) 


(R1-M1)(R1-M1)T 
e  j  (r2-m2)(r1-h1)t 


(R1-H1)(R2-H2)T  ...  (R-j-M^ 

(R2-M2)(R2-M2)T  ...  (R2-M2)(VMm)T 

•  • 

•  • 

•  • 

(R.-VO'zV  -  <W‘WT 


where 


It  is  apparent  that  m^.  =  ai£  under  hypothesis  Ha  and  m-j£  *  b^£ 
under  hypothesis  Hk.  For  an  m-element  array,  the  mz  submatrices 
may  be  classified  into  two  types:  the  one  associated  with  the  ith 
channel,  namely,  the  autocovariance  submatrix 

(2-22)  A  E  {(R.  -  Mi)(RrM1)T} 

and  the  crossvariance  submatrix 

(2-23)  A  E{(R1-M1)(Rj-Hj)T>  ,  1  f*  j. 
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which  Is  associated  with  the  1th  and  jth  channels,  respectively. 
These  submatrices  may  also  be  expressed  in  terms  of  the  elements 
of  Rj  in  Eq.  (2-15),  as  follows, 

(2-24) 


VE{ 


(rn-n,ii)(rji-mji) 

|(  r1 2”m1 2 )  ^  rj  1  “m<j  1 ) 


( rfl“mi  1 )  ( rj2“mj2 ) 
( ri  2“mi  2 )  ^  rj  2"mj2  ^ 


(ril-mil)(rjN-mjN) 
( ri  2”m1 2 )  ^  rj  N"m  j  N ) 


(f1N 


-m 


IN 


>(rjrV 


^iN^iN 


^rj2"mj2) 


riN"miN 


We  proceed  to  examine  the  elements  of  Kjj  under  both  hypotheses. 
Since  rfe  and  rjn  are  statistically  independent  Gaussian  random 
variables,  and  processes  n.(t)  and  w. (t)  are  uncorrelated,  we  have 
under  H  1  1 

O 


(2-25)  E{(ri£-mu)(rjn-mjn)|Ha)} 

E  { (ri£rjn)|  HaJ  -  a^3jn 

a  |J  T  iCt-u+U-IJv]  •  uithnM  .  dudt  +  5£n*  hw 

-  V»n  •  \  *  <Vh„>  •  j.T  ♦,<*>  •  ♦.fHU-D.D  •  dt 

The  elements  of  for  1j*j  hence  become 

fT 

(2-26)  EUru-mu)(r;jn-mjn)|Ha}  =  (xn-hw)  •  j  ^(tW^t+U-iJ^dt 


and  the  elements  of  Kti  become 


(2-27)  E((ru-mu)(r1n-inin)|Ha} 


=  ‘W  '  j.T  **(t)  +n(t)dt  +  6< 


if  a  f  n 
if  a  =  n 


The  elements  of  K-ji  and  Ki i  under  hypothesis  Hb  may  be  derived 
in  the  same  manner.  It  is  found  that  they  are  identical  to  those 
under  hypothesis  Ha.  Consequently,  the  covariance  matrices  under 
the  hypotheses  defined  in  Eqs.  (2-7)  and  (2-8)  are  equal,  and  the 
hypothesis  detection  must  rely  on  the  difference  in  the  mean  M 
associated  with  each  hypothesis.  The  autocovariance  submatrix 
and  the  cross  covariance  submatrix  have  the  following  forms: 

fx,  0  0  0  *1 


(2-28)  Ki.  = 


0  X- 


(2-29) 


0  0  0  X, 


(»rV  }  T  ♦j(t)*|Et*(J-1)v]dt  (yh*)  ♦1(t)t2Ct+(J-1)v3dt  •••  (yhj  |TT«1(t)eNrt+fj-i)v]dt 

(W  |t  (yy  ytjyXJ-OvDdt  •  ••  (yhj  j\2(t),Nrt+(j-i)v]dt 

•  *  . 

•  •  • 

■  •  • 

_(»rhw)  |  t  ♦N(t)»1Ct+(J-1)v3dt  (x2-hw)  ...  (xN-hw)  |\N(t)<N[t*(j-i)v]dt 


The  element  of  Kjj  at  the  an  position, is  equal  to  that  of  K.. 
at  the  na  position,  because  J1 


11 


% 


(2-30)  (A^hJ  j  *A(t)  .  *nCt  +(j-1)v:  •  dt 

■  |  ^(u)  •  J  RCt+( j-1  )v-u3  •  *A(t)  .  dt  du 

=  (xi  "  V  |  *  ♦4CU+(1“J)V3  •  du* 

Hence,  Kij=Kl-j  and  consequently  the  covariance  matrix  K  is  symmetric. 
In  terms  of  K^.  and  K^,  the  covariance  matrix  K  may  be  written 

K11  K12  Klm 

(2-31)  K  =  ^  K22  ... 


“ml  ^2  ***  Sim 


Since 

(2-32)  Q  A  K"1  =  (K1)"1  -  QT 
the  Q  matrix  Is  also  symmetric. 

C.  Broadside  Interference 

To  compute  the  Q  matrix,  one  has  to  go  through  a  tedious  process 
First,  the  Integral  equation  defined  by  Eq.  (2-11)  has  to  be  solved. 
Afterwards,  one  has  to  compute  many  Integrals  In  order  to  get  K. .. 
Finally,  one  must  Invert  the  covariance  matrix  K.  1J 

Referring  to  Eq.  (2-29),  we  see  that  a  much  simpler  expression 
for  K  Is  obtained  If  the  external  Interference  arrives  from  the 
broadside  direction  of  the  array,  l.e..  If  v  *  0.  When  this  Is  the 
case,  the  orthogonality  property  of  {♦£(t)>  ensures  that  all  off- 
diagonal  elements  of  K^h,  1j*j,  vanish.  Furthermore,  all  cross¬ 
covariance  submatrices  Become  equal,  as  follows: 


0 


0 


(2-33)  j(v=0)  = 


Arhw 

0  x2"hw 


0  0 


Vhw 


±Kh 


Let  K*  be  the  covariance  matrix  for  this  special  case.  Then  the 
expression  for  each  element  of  Q*,  defined  by  Q*  =  K*-1,  can  be 
obtained  by  generalizing  the  results  for  the  inverses  of  the 
covariance  matrices  associated  with  a  two-element  array,  a  three- 
element  array,  etc.  It  is  found  that  all  submatrices  Q*.  are 
identical  and  diagonal.  The  nonvanishing  elements  11 


<*«’  %+0(N+£)*  2{2N+£)(2Nn)*  etc* 
may  be  written  as+ 


+  (m-2)(x£-hw) 


<2‘34)  ,f(1-1)N+tX(,-l)N+e]'x^)xt  •  (xrhw)  -  -hj* 


'a  ’V 


A£  +  (m-2)  .  (x£-hw) 
hw  *  "  (m-1)hwD  * 


1  <  i  <  m 

1  <  l  <  N 


The  submatrices  Qfj,  1j*j,  are  equal  and  diagonal  as  well.  The 
diagonal  elements  may  be  expressed  as 


<2'35)  "f'1-'  )N+o  •  m 

Combining  Eqs.  (2-34)  and  (2-35),  we  have 
<2'36)  ^ 


+A  simple  proof  of  Eq.  (2-36)  Is  shown  in  Appendix  A.  Notice  that 
Q*j  t  K*j '  for  all  1  and  j. 


Although  Eq.  (2-36)  is  only  good  for  v  =  0,  the  result  of  this  special 
case  can  be  generalized. 

D.  Interference  from  Arbitrary  Angles 

Although  the  digital  computer  may  be  used  to  obtain  the  Q  matrix 
numerically  If  the  order  of  K  is  not  too  large,  we  hope  to  obtain  a 
general  analytical  result.  Thus  we  shall  have  to  find  another  method 
for  computing  the  Q  matrix. 

The  simplicity  for  the  case  of  broadside  interference  is  due  to 
the  Interfering  signal  having  equal  phase  at  each  element.  Suppose 
we  consider  the  following  delayed  signal 

(2-37)  r*(t)  a  ta.(z)  +  n^z)  +  w.(z)}z=t_jm_i^ 

«  aCt-(m-i)P-(i-l)xIl  +  nCt-(m-l)v]  +  w.Ct-(m-i)vH 

where  Eq.  (2-3)  has  been  used  to  replace  ai(z)  and  n-j(z).  The  second 
term  in  Eq.  (2-37)  does  not  have  the  index  i,  becuase  the  interference 
has  the  same  phase  for  all  r|(t)'s.  We  examine  the  covariance  matrix 
associated  with  r*(t).  Let 

(2-38)  r*£  A  |  r|(t)  *£(t)  •  dt 


(2-39)  m*£  a  E{rt  |HR>  . 


Then  the  element  of  the  Ijth  submatrix  at  the  &n  position  is 
(2-40)  E((rTl-^t)(rJ„-*n)) 

*  ||  E{n[t-(m-l)vD  •  nCu-(m-l)v]}  •♦£(t)«|>n{u)  •  dt  du 
+  |J  ElWjCt-On-iM  •  WjCu-(m-j)v]}  •  4ft(t)4n(u)  *  dt  du 


=  Whw>  + 


6Ct-u+(i-j)v]  •  ♦£(t)«n(u)  .  du  dt 


-  *< 


6zn  *  ^xn  "  ^ 


azn  *  xn 


if  i-j 


14 


Equation  (2-40)  indicates  that  the  covariance  matrix  of  r^(t)  is 
identical  to  that  for  broadside  interference.  This  is  also  true  if 
a-j(t)  in  Eq.  (2-37)  is  replaced  by  b. (t),  which  is  the  signal  under 
hypothesis  Hfa.  1 

* 

Physically,  r,*(t)  is  the  signal  obtained  by  delaying  the  signal 
r-j(t)  in  the  ith  channel  by  t-(m-i)v  seconds.  In  general,  such  time 
shifting  of  the  received  signal  could  be  done  with  a  delay  line 
behind  each  array  element.  Since  the  w^(t)  are  assumed  to  be  wide- 
sense  stationary  processes,  shifting  on  the  time  axis  will  not 
affect  their  statistics.  We  next  examine  how  the  rt.  's  are  related 
to  the  r.£'s.  We  know  that  u 

N 

(2-41)  r.(t)  =  lim  £  r.  <{>.(t),  -T  <  t  <  T 

1  N-*°  .1=1 

It  is  also  true  that 

N 

(2-42)  r.Ct-(m-i>]  =  lim  l  r*  Tt-(m-i)v],  -T  <  t-(m-i)v  <  T 

1  N-~  £=1  1£  £ 

and 

N 

(2-43)  rt(t)  A  r.[t-(m-i)v3  =  lim  £  r*<*(t) 

1  1  £=1  1)1  1 


Equating  Eqs.  (2-42)  and  (2-43),  we  have 

(2'44)  j,  =  j,  r1« 

If  both  sides  of  Eq.  (2-44)  are  multiplied  by  <j>n(t)  and  integrated 
between  the  limits  -T  and  T,  the  result  yields 

(2-45)  r|n  =  l  ri£  •  j  $n(t)  •♦jJXni-Ov]  •  dt 
£—  1  — T 


We  define  the  vector  R* 
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i 


where  0  is  a  null  matrix  and  the  submatrix  C.  is  a  square,  un- 
symmetric  matrix,  defined  by 

(2-49) 


a 


•  *2Cv-(m-1)v]dt  ••• 

Cl  * 

J  *2(t)  • 

j:t  *2^) 

•  ♦2[t-{n-i)vDdt  ••• 

J  <2(t)*<‘N[t-(m-i  )vZ<Jt 

(  $N(t)  •  +j[t-(m-i)v3dt 

iv> 

*  ♦2[t-(m-i)v]dt  ••• 

fT 

J  eN(t)*$Nrt-(n-i)v]dt 

In  matrix  notation,  Eq.  (2-45)  may  be  expressed  as 


(2-50)  lim  R*  =  lim  CR 
N-x»  N-*» 

We  have  pointed  out  earlier  that  if  we  compute  E{(R*-M*)(R*-M*)T}, 
where  M*  a_  E{R*|H),  we  get  the  same  covariance  matrix  as  if  the 
interf erring  source  were  at  broadside.  In  other  words. 


(2-51)  K*  =  E((R*-M*)(R*-M*)  > 

Substituting  Eq.  (2-50)  into  Eq.  (2-51),  yields, 

(2-52)  K*  =  CKCT 

From  Eq.  (2-52),  it  can  be  shown  that 


(2-53)  Q  =  C1  Q*C 

Eq.  (2-53)  provides  us  with  a  much  easier  way  to  compute  the  Q 
matrix.  Since  Eq.  (2-42)  is  defined  for  -T  <_  t  -  (m-i)v  <  T  and 
the  series  given  by  Eqs.  (2-41)  and  (2-43)  converge  to  rift)  and 
rf(t)  respectively  if  N  approaches  infinity,  Eq.  (2-53)  is  correct 
and  exact  if  those  two  conditions  hold.  For  a  finite  N,  Eq.  (2-53) 
is  merely  an  approximation. 
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1 


Now  we  use  this  transformation  to  find  the  elements  of  Q  for 
an  arbitrary  <p.  Performing  the  matrix  multiplications  indicated 
by  Eq.  (2-53),  we  find  that  the  submatrix  Q..,  which  is  one  of 
m2  of  its  kind,  is  given  by  J 

(2-54)  Qij  =  Q*j  cj» 

If  the  last  term  for  the  series  expansion  is  denoted  as  N,  then 
there  are  N2  elements  in  the  square  matrix  Q-ji.  We  use  £,  with 
proper  subscripts,  to  denote  elements  of  Q.  From  Eq.  (2-54),  we 

find  the  element  of  Q..  at  the  £n  position  is  equal  to 

i  J 

(2-55)  ^ 

aC(1-l)N+tX(j-l)N+nD  =  •  ♦nCu-(i«-J)v> 


n 

‘  (1-1  )N+kX(j-l )N+k]  *  <t>k(t)<J>k(u)'dtdu 


Substituting  Eq.  (2-36)  into  Eq.  (2-55),  yields 


(2-56) 


«ii  fT 

3C(i-l)N+£X(j-l)N+n]  =  J-T  *£nt-(m-i)v><frnCt-(m-j)v]  •  dt 

“X  ^Y-fm-1  *  |j_T  ^Ct-(m-i)vD  .  4»nCu-(m-j)v> 


*  4>k(t)  <J>k(u)  •  dt  du 


We  have  made  use  of 


(2-57)  6(t-u)  =  I  <|>,(t)  ♦.(u) 

£=1  *  * 

to  get  Eq.  (2-56).  This  result  will  be  used  to  formulate  the 
analytical  expression  for  the  test  statistics  in  the  next  chapter. 
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CHAPTER  III 

THE  OPTIMUM  DETECTOR  FOR  DETECTING  COMPLETELY 
KNOWN  BINARY  SIGNALS 


A.  The  Likelihood  Ratio  Test 


In  this  chapter,  the  equations  of  the  test  statistic  and  those 
required  for  showing  the  receiver  operating  characteristic  of  the 
optimum  processor  are  formulated.  The  results  presented  here 
optimize  the  performance  of  the  entire  detection  system,  rather 
than  the  antenna  array  or  the  receiver  separately. 

The  binary  signals  a(t)  and  b(t)  are  assumed  known  exactly. 

Thus  the  uncertainty  in  the  signal  at  the  input  of  each  channel 
arises  solely  from  the  additive  random  noise  which  is  the  sum  of 
the  interference  and  the  internal  noises.  The  detection  system 
will  have  to  make  a  choice  between  two  hypotheses  by  processing 
the  signal  from  each  channel.  The  signal  to  be  processed  at  the 
ith  channel  can  have  one  of  two  different  signals: 

Hypothesis  H,:  r,(t)  =  a-(t)  +  n.(t)  +  w.(t) 

a  I  I  i  1 

Hypothesis  Hfa:  r^(t)  =  b,.(t)  +  n.(t)  +  w.(t) 
where  1  £  i  $  m. 

Dual  hypotheses  testing  relies  on  a  decision  rule  for  dividing 
the  observation  space  into  two  parts.  Various  optimum  decision  rules 
(e.g.,  Bayes  criterion,  Neyman-Pearson  criterion,  etc.)  can  be  used. 
These  decision  criteria  for  the  optimal  receiver  all  lead  to  a  like¬ 
lihood  ratio  test  (21j: 

p(R|Hb)  Hb 

(3-1)  L.R.T.  A  lim  u  \  >  threshold 

N-**»  PvK  < 


where  N  is  the  last  term  of  the  Karhunen-Loeve  expansion.  The 
structure  of  t!e  optimum  detector  immediately  follows  once  the 
result  for  the  likelihood  ratio  test  is  obtained. 

The  probability  density  function  of  R  under  hypothesis  H 
for  a  N-term  Karhunen-Loeve  expansion  and  an  m-element  array  ts 
equal  to 
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1 


i 


% 

% 

> 


i 


i 


V' 


A 


-  i  (R-A)TQ(R-A) 

(3-2)  p(R|Ha)  =  - e  2 

^r)mN|K| 

The  exponent  may  be  written  as 

(3-3)  -  \  (R-A)T  .  Q  •  (R-A) 

1  m  m  N  N 

=  “7^  j-I-,  nl}  {rW&u)  ’  -aCd-DN+OCU-DN+n]*  (rjn“ajn) 

“  -  7  Jl,  |^TCr1(t)-«i(t)Xrj(2)-aj(z): 

N 

*  1 1  £=1  3C(i-l)N+£X(j-DN+n3  *n(z)  '  dt  dz* 

We  may  express  Eq.  (3-3)  in  a  more  compact  form  by  defining 

(3  4)  Hij(t,z)  Lll^  +n(z) 

Since  the  Q  matrix  is  symmetric,  it  can  be  shown  that  H..(t,z)  = 
H^(z,t)  and  thus  J 

m  rrT  n 

(3-5)  I  I  jj  D^OOa^z)  +  aj(t)  rj(z)]  .  Hjj(t,z)  •  dt  dz 
i , j —I  — T 

'  2  ^  .^  {}  T  r -j (T)aj(z)  •  H^t.z)  •  dt  dz 

u 

Hence,  in  terms  of  H^t.z),  the  exponent  becomes 


(3-6) 


-  \  (R-A)T  •  Q  •  (R-A) 

=  -i  ll  ff  r-(t)r.(z)  *  H^.(t,z).dtdz 
c  i ,3=1  JJ-T  J  13 

■i  [  I  [f  a.(t)  a.(z)  •  H? .(t,z)  •  dt  dz 

c  i ,j=l  JJ-T  1  3  13 

+  I  I  ff  r.(t)  a,(z)  •  H^-(t,z)  •  dt  dz 
i,j=l  JJ-T  1  3  13 


The  probability  density  of  R  under  is  readily  obtained  by 
changing  a(t)  in  Eq.  (3-6)  to  b(t).  We  compute  the  likelihood 
ratio 


(3-7) 


L.R. 


lim 

N-*» 


p(R|Hb) 

p(R|Ha) 


=  lim 
N-*» 


l  ^  fj  T  ^(t)  •Cbj(z)-aJ.(z)]«  H^(t,z).dt  dz 

m  f  r T  M 

II  [b.(t)b.(z)-a.(t)a.(z)>  H^,(t,z).dtdz  f 

i,j=1  JJ-T  ‘  3  13  13 


ri(t),  the  input  of  the  ith  channel,  appears  only  in  the  first  term 
of  Eq.  (3-7).  Since  a(t)  and  b(t)  are  assumed  known  exactly,  the 
second  term  thus  is  a  constant  and  may  be  incorporated  in  the 
threshold.  Let 


(3-8) 

Then, 

(3-9) 


s(t)  =  b(t)  -  a(t) 

from  Eq.  (3-7),  an  equivalent  likelihood  ratio 


m 


ELRT  A  J  I  ri(t)si(z)  •  H~  (t,z)  . 

t,j=l  JJ-T  1  3  13 


test  (ELRT)  is 
dt  dz 


Sb 

<  threshold 


where  we  have  taken  the  limit  as  N+*>: 


(3-10) 


H^Ct.z)  =  lim  Hj.(t.z) 

N->oo  'J 

N 


=  N~  Lnll  ^O-DN+Oaj-DN+n]  *i(t)  *n(z) 
fT 

=  TT^"  <sCt-t,+(m-i)vH  •  sCz-t'+On-jJvI]  •  dt1 
nw  J-T 

-  k  j,  lfVCt-t'+(m-Dv]  ’ 


•  6Cz-u'+(m-j)vD  •  ^k(t ' )  ^(u1)  •  dt1  du1 
Substituting  Eq.  (3-10)  Into  Eq.  (3-9),  yields 


ELRT 

i  m  rT 

=fiw  ill  J.T 

where 

m  r  rj 

-  n 

1,j=l 

rri<V5j<zj>  • 

(3-11) 

tj  A.  t-(m-1  )v 

(3-12) 

z j  A  z  -  (m-j)v 

and 

(3-13) 

,  %  •  *AthAz) 

y(t,z)  A  l  -i - 1 - 

*=1  _  .  "w 

m  +  t - L- — 

i  w 


Let  n  be  the  threshold.  The  optimum  detection  system  Is  defined  by 


(3-14) 


m 

L(R)  A  l 
i=l 


m 

n 

i  »j=l 


fT 

Trl(tl)sj(*j)  *  1'(t»z)  *  dt  dz 


Sb 


*1 


H 


a 


The  likelihood  ratio  test  leads  us  to  the  test  statistic 
L( R) ,  which  can  be  generated  by  correlating  the  delayed  signal 
r-f(t-f)  with  some  known  functions.  We  compare  L(R)  to  a  threshold 
to  make  the  decision. 


In  Eg.  (3-14),  the  known  signal  s(t)  and  the  array  received 
signal  r.(t)  are  time-shifted,  because 


(3-15)  Sjftj)  -  s.Ct-(m-i)vU 

=  b[t-(i-l)T- (m-i)vl  -  a[t-(i-l)T-(m-i)v] 
and 

(3-16)  ^.(t.)  =  r.[t-(m-i)v3 


The  delay  time  (i-l)x  occurs  because  each  array  element  has  unequal 
distance  to  the  signal  source.  But  the  additional  delay  (m-i)v 
results  from  computing  the  Q  matrix  using  the  transformation  C. 

The  double  summation  term  in  the  test  statistic  L(R)  is  influenced 
by  the  power  spectrum  of  the  interference,  for  <{,(t,z)  is  obtained 
from  the  integral  equation  defined  by  Eq.  (2-11).  The  double 
sumnation  term  may  be  rewritten  as 


(3-17) 


where 


(3-18) 


s^)h(z)  •  dz 
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Hence,  is,  in  general,  correlated  with  m  distorted  known 

signals  si(zj)  because  sj  's  are  weighted  unequally.  For  the  case 
when  the  interference  isJa  white  noise  process  with  autocorrelation 
function  R(t-u)  =  hn  .  6( t-u) .  it  is  found  by  solving  Eq.  (2-11) 
that  A^hj^  —  hn  for  all  £.  Thus  the  weights  for  the  st^'s  are 
identical  and  the  test  statistic  takes  a  much  simpler  form: 


(3-19) 

m  fT 

L(R)  s  I  r.(t.)s.(t.)dt  - 
i=l  J-T  1111 


I  I  ri(ti)si(t1)dt 
i ,3=1  J-T  1  1  J  3 


A  scheme  for  processing  the  input  of  the  ith  channel  based 
on  Eq.  (3-14)  is  shown  in  Fig.  3-1. 


Fig.  3-1— The  ith  channel  processing  unit 


B.  Receiver  Operating  Characteristic 

Since  the  test  statistic  L(R)  is  obtained  by  the  linear 
operation  indicated  by  Eq.  (3-14)  on  rj(tj),  it  is  a  Gaussian 
random  variable  under  either  hypothesis.  The  density  of  L(R)  is 
found  once  its  mean  and  variance  are  determined. 

The  mean  of  L(R)  is  readily  obtained  if  we  replace  r.(t. )  in 
Eq.  (3-14)  by  its  expectation,  i.e., 

(3-20) 

mLa  a  E{L(R)  |Ha) 

m  fT  m  rT 

=  l  a.UOs^tOdt  -  l  l  ai(ti)s1(zi)*y(t,z)*dt  dz 
i*l  J-T  1  1  1  1  i,j=lJ'-T  1  1  J  J 

and 

(3-21) 

mLb  a  E{L(R)|Hb} 

m  ,T  m  r  ,T 

*  I  MtJs.ftJdt  -  [[  b’ftjs.fzj-^t.zj.dt  dz 
i=l  i-T  1  1  1  1  i td=l  1  1  J  3 

The  variance  of  L(R)  under  either  hypothesis  may  be  calculated  from 
the  following  relation 

(3-22)  var{L(R)|Hk}  =  E{L2(R)|Hk>  -  E2{L(R) |Hk> 


Evaluating  the  variance  is  straightforward  but  tedious.  The  derivation 
is  given  in  Appendix  B.  We  find  that  the  variance  of  L(R)  equals 
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By  expressing  s^t^)  as, 

S1  <ti )  = 

2 

It  can  be  shown  that  the  variance  of  may  be  written  in  terms  of 
the  distance  between  the  means  of  tne  test  statistic  under 
either  hypotheses  as 

(3-24)  o*  -  hw  (mLb  -  m^) 

2 

Since  of  is  nonnegative,  one  concludes  that 
(3-25)  mLb  >  mLl 

We  next  formulate  the  equations  for  the  probabilities  of  false 
and  correct  detection.  Let  the  probability  of  choosing  %  when  Ha 
is  present  be  denoted  by  Pc  and  that  of  choosing  %  when  %  is 
present  by  Pq.  Then,  as  illustrated  in  Fig.  3-2,  the  probability 
of  false  detection  1s+ 

(3-26)  Pp  =  J  p(L(R)|Ha)  •  dL(R) 


and  the  probability  of  correct  detection  is 


+erfc*(x)  is  defined  as 


dy 
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(3-27) 


,CO 


p 


0 


p(L(R)|Hb)  •  d  L(R) 

(L-mLb)2 


1 

127 


e 


dL  =  erfc* 


The  Pp-Pn  plot  1S  referred  to  as  the  receiver  operating  characteristic 
(ROC). 


Fig.  3-2— Probabi  1  i  ty  density  functions  of  L(R). 


For  a  communication  system,  we  are  usually  interested  in  the 
total  probability  of  making  incorrect  decisions.  If  the  a  priori 
probabilities  of  the  source  are  equal,  i.e.,  if  PrCHaJ  =  PrQib3  = 
then,  Pe,  the  total  probability  of  making  error  is  simply 

(3-28)  P£  =  %  PF  +  %  (1-PD) 

For  a  given  mLa,  mih  and  o(_,  the  probability  of  error  varies  with 
the  threshold  n.  The  threshold  n0  which  minimizes  Eq.  (3-28)  may 
be  determined  by  solving 


dP 

(3-29) 

II 

o 

for  n. 

The  solution  is 

(3-30) 

no  =  4  (mLb  +  "’La5 
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The  corresponding  probability  of  error,  P  (n. ),  Is  found  to  be  equal 
to  f  0 


(3-31)  Pe(n„)  -  e,fc,(yM  ) 

which  Implies  that  Pe(nn)  may  be  minimized  by  maximizing 
(mi.b-niLa)/0L«  Alternatively,  by  substituting  Eg.  (3-24)  into 
Eq.  (3-31),  we  obtain  another  expression  for  P  (tv). 


(3-32)  Pe(n0)  *  erfc*^ 


mLb-mLa 


1 

1 


CHAPTER  IV 

THE  DIRECTIONAL  CHARACTERISTICS 
AND  THE  ERROR  RATE  PERFORMANCE 


A.  Introduction 


To  make  the  general  theory  developed  in  the  last  chapter  clearer, 
we  shall  treat  the  case  when  the  binary  signals  are  of  the  form. 


(4-1)  Hfl:  a(t)  =  ea(t)  •  cos(ut  +$  &) 

(4-2)  Hb:  b(t)  =  eb(t)  •  cos(oit  +4,  $) 

The  signals  given  above  represent  the  amplitude-shift  keying  (ASK) 
binary  signals  if  the  modulating  signals  ea(t)  and  e[j(t)  only 
differ  in  amplitudes.  If  the  modulating  signals  are  related  by 
ea(t)  =  -  eb(t),  then  the  signals  given  in  Eq.  (4-1)  and  Eq.  (4-2) 
are  biphase  modulated  binary  signals. 

The  optimum  communication  system  contains  an  array  and  an  optimum 
signal  processor.  The  task  for  the  system  is  to  detect  the  binary 
signals  arriving  from  angle  e  in  the  presence  of  an  interfering 
signal  from  the  given  angle  $.  To  accomplish  this  goal,  we  in¬ 
tuitively  expect  the  system  to  function  in  favor  of  the  signal 
arriving  from  the  desired  angle  e.  To  investigate  this  aspect,  we 
study  how  the  system  treats  a  testing  signal  arriving  from  an  arbitrary 
angle  e.  Since  the  detection  system  is  based  on  a  comparison  of  the 
processor  output  L(R)  with  a  threshold,  we  examine  L(R)  to  determine 
the  system  characteristics.  One  important  characteristic  of  L(R)  is 
its  mean.  Hence,  we  study  the  mean  instead.  The  performance  of  a 
detector  is  judged  mainly  by  its  ability  to  make  correct  decisions. 

For  this  reason,  we  shall  also  discuss  the  error  rate  performance  and 
the  methods  of  optimizing  it. 

Here,  we  shall  assume  that  the  waveform  for  the  testing  signal 
is  the  same  as  the  signal  a(t).  The  reasons  for  this  choice  are 
(1)  to  simplify  the  mathematics  involved,  and  (2)  to  understand  the 
behavior  of  the  mean  mi_a  through  studying  the  directional  character¬ 
istics  of  the  optimum  system.  The  output  from  the  optimum  detector 
( Q,4>)  when  a  signal  a(t)  arrives  from  angle  e  is  given  by* 


ft 

We  use  "detector  ( e ,4>) "  to  mean  a  detector  designed  to  detect 
the  desired  communication  signal  arriving  from  angle  0  in  the 
presence  of  the  interfering  signal  with  incident  angle  $. 


where  ^(t^)  Is  defined  as 

^(tj)  in^^)  +  wi(ti)  +  ai(t1) 

a  n^tj)  +  w^t.)  +  aCt  -  (i-1)  £  cos  e  -(m-i)v] 

The  mean  of  the  system  output  L(R)  is 


(4-4)  mu(e,*;e)  A  E{L(R)> 

in  rT  A 

=  I  a.D>(m-i hi  •  s.Ct-(m-i)v]  •  dt 
1=1  L  t  1  1 

-  I  l  [f  a.[t-(m-i)v]  •  s.Cz-(m-j)v]  • 

1,j-l  jj-T  1  3 

.  4'(t,z)  •  dt  dz 

We  write  mLa  as  a  function  of  the  angles  to  emphasize  its  dependence 
on  these  angles.  The  m.  appearing  in  the  last  chapter  is  just 

mLa(e.*;e). 

The  closed  form  representation  for  mia(e,<|>;e)  will  be  derived  in 
the  next  section.  In  Sections  IV.C  and  IV. E,  the  directional  char¬ 
acteristics  and  the  error  rate  performance  are  investigated.  The 
discussions  are  restricted  to  the  case  when  the  interfering  signal 
is  a  white  noise  process  and  the  desired  signals  are  biphase  modu¬ 
lated.  Section  IV. F  generalizes  these  results  to  the  case  when  the 
Interference  is  a  general  Gaussian,  wide-sense  stationary  process 
with  arbitrary  power  spectrum. 


B.  The  Mean  of  L(R) 


We  assume  that  the  frequencies  contained  in  the  modulating  signals 
ea(t)  and  e^t)  are  much  lower  than  the  carrier  frequency  so  the 
following  approximation  may  be  made* 


(4-5)  eki Ct-(m-1  )p]  £  ek(t),  k  =  a  and  b. 


Hence,  a..Ct-(m-i )v]  may  be  approximated  by 


(4-6)  a.Ct-(m-i)v]  %  ea(t)  •  cos{u>[t-(m-l  )v]  +  *  -(i-1  )w(x-v)  }. 
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In  the  same  manner. 


(4-7)  Ct- (m-i ) vD  %  e(t)  •  cos{a£t-(m-l  )v]  +  -( i-1  )w(t-v) }, 


where  e(t)  A_  eb(t)  -  ea(t). 
By  making  use  of  the  identities 


sin(  x)  •  cos(  x 


III  3  I  ll  I  *  A  I  • 

l  cos(i-l  )x  =  — i 
1=1  sin(|) 


m  sin/2-x)»  sin(^^x) 

I  sin(i-l)x - L If - IX. ± 

i=l  sinf£j 

it  can  be  shown  from  Eq.  (4-6)  and  Eq.  (4-7)  that 


♦This  approximation  introduces  a  maximum  time  shift  of  d/c  (m-1) 
seconds  in  the  modulation  envelope.  We  assume  this  time  shift 
to  be  negligible  at  the  maximum  modulation  frequency. 


i 


(4-8)  l  s.[t-(m-i)v] 
i=1  1 


sin  ■j-  w(t-v) 
sin  ^  (t-v) 


•  e(t)  •  cos{u>Ct  -  — 2”  (t+v)3+ 


.  j 

i 

j 


m  A 

(4-9)  l  a.Ct-(m-i)v] 
i=l  1 


and 


(4-10) 


m 


sin  V|-  to(x-v)  ,  A 

- — - - -  e  (t)*cos{a>[t  — k—  (x+v)3  +  } 

sin  £  (?-v)  a  2  s 


m 


^  a.Ct-(m-i)v]  •  s.Ct-(m-i)v3 


e(t)e  (t)  sin  £  w(t-t)  •  cos  ^r-  u(t-t) 
^  sin  £  (t-t) 

e(t)ea(t)  sin  £  w(t-t-2v) 


m-1 


sin  £  (t+t-2v) 


cos{2w(t  -  ^nv)+2(j)g  -  (t-t)} 


Substituting  Eqs.  (4-8),  (4-9)  and  (4-10)  into  (4-4),  we  obtain 

(4-H)  _  A 

A  sin  y  oj(t+t-2v)  fT  e  (t)e(t) 

- - - - 4 - 

sin  £  (t+t-2v)  '  -T  * 


mLa(e,<j>;e)  = 


•  cosC2w(t  -  — 2^*  v)  +  -  ~2~  w(t+t)D  •  dt 


(4-11) 

(Cont.) 


sin  ijj(t-t)  •  cos  **»(T"T) 

■T  ea(t)  e(t) 

•  to  /  A  \ 
sin  -p  (t-t) 

-T  2 

sin  SJ-  w(t-v)  sin  Sy  w(t-v) 

X  < 

c  ^  .  y  i 

sin  j  (t-v)  sin  £  (t-v)  £.=1  m  ,  ~ 

\  -hw 


ea(t)e(z)*‘ti(t)$Jt(z)*cos{uCt  -  — (t+v)]  +  <>sl* 


•  COS{u£z 


(t+v)3  +  4>s } 


dt  dz  . 


If  uis  large,  the  first  term  in  Eq.  (4-11)  is  comparatively  small 
and  may  be  neglected  for  computation. 

When  the  interference  is  a  Gaussian  white  noise  process  with 
spectral  height  hn,  Eq.  (4-11)  reduces  to  a  very  simple  form. 


(4-12) 


mLa(e,4>;e)  * 


T  e(t)ea(t>  dt  rn_  m-1  , 

- 2 -  dt  »  COS  — ij—  u>( t-t) 


sin  S^  w(t-t)  ^ 


sin  Tj-  w(x-v)  sin  S^  oj(t-v) 


sin  ip  (t-t)  m  +  "jjw  S^n  ^  sin  ^  (t-v) 


where  numerically  small  terms  in  Eq.  (4-11)  have  been  neglected. 
C.  The  Directional  Characteristics 


The  detection  system  is  assumed  to  operate  with  two  signal 
sources,  the  desired  communication  signal  and  the  interference. 

The  incident  angles  for  these  sionals^are  e  arid  <t>,  respectively. 
Hence,  the  behavior  of  mLa(e,4i;e)  at  0=0  and  e=<j>  interest  us  most. 

For  this  reason,  we  next  study  mLa(0»4>*»0)  and  mLa(e,<K4>).  In  this 
section,  the  interfering  signal  is  assumed  to  be  a  Gaussian  white 
noise  process  of  spectral  height  hn< 

The  test  function  L(R)  is  obtained  by  processing  the  sum  of  the 
desired  signal,  the  interference  and  the  internal  noise  in  the  manner 
given  in  Eq.  (3-14).  The  quantity  mL  (e,<j>;o)  is  the  mean  of  the 
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processor  output  L(R)  under  hypothesis  Ha.  Since  the  interference  and 
the  internal  noise  have  been  assumed  to  be  zero-mean,  mLafe.^e)  in 
fact  is  equal  to  the  component  contributed  to  L(R)  by  the  signal  a(t) 
from  the  desired  sQurce  for  wjiich  the  processor  is  designed  to  be 
optimum.  Letting  e*e  (i.e.,  t=t)  in  Eq.  (4-12),  we  obtain 


(4-13)  mLa(e,4>;e) 


1 


2  m  I  x  ' 

sin  Tf-  oi(t-v) 


_  .  "w  sin 
m+^ 


(x-v) 


Hr 


ea(t)e(t) 


dt 


For  the  same  reason,  mj_a ( © •  4> i 4* )  is  equal  to  the  component  contributed 
to  L($)  by  the  signal  a(t)  arriving  from  the  angle  identical  to  that 
of  the  interference.  Letting  6=<j>  (i.e.,  t=v)  in  Eq.  (4-12),  we  find 
that 


(4-14) 


mLa(0*  =  cos  TT” 


•T  ea(t)e(t) 


It 


dt 


sin  y  t»»(r*v) 
sin  y  (t-v) 


1  - 


m 


m  + 


w 

* 


It  can  be  shown  that  mLa(e,<t>;e)  never  changes  sign  for  a  given 
fT 

j  ^  ea(t)e(t)  .  dt 

as  the  angles  of  arrival  or  the  number  of  elements  are  varied.  More 
precisely,  the  sign  of  mLa(e,<|>;e)  is  always  the  same  as  the  sign  of 

fT 

j  T  dt* 

It  can  be  shown  from  Eq.  (4-14)  that  mLa(e»<f>;<f>)  vanishes  if* 


Equation  (4-15)  is  derived  from  setting  sin  £  a>(T-v)a0.  x  Is 
the  free-space  wavelength  associated  with  the  carrier.  The 
quantity  m,  *(e,4>;4>)  also  vanishes  if  cos  !2ii  w(t-v)  =  0. 

La  2 
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(4-17) 

0  1  |i\a(0.<f>;4>)  I  i 


Since  the  upper  bound  for  Jmi  a(e*4>;4>)  I  is  just  the  lower  bound  for 
|mLa(0»4>;0)  | »  the  former  will  never  become  greater  than  the  latter. 
They  are  equal  only  when  0=<t>. 

The  bounds  given  in  Eq.  (4-16)  and  Eq.  (4-17)  indicate  that  the 
processor  produces  unequal  outputs  for  the  signal  a(t)  when  this 
signal  arrives  from  different  angles.  More  precisely,  the  processor 
functions  in  favor  of  the  signal  a(t)  when  it  arrives  from  the 
desired  angle  0  and  suppresses  it  if  the  incident  angle  is  <p. 

For  a  fixed  m  and  d,  the  ability  of  each  processor  to  suppress 
the  signal  a(t)  arriving  from  the  angle  <f>  depends  on  the  angular 
separation  between  the  incident  angles  e  and  <f>.  To  show  this,  two 
typical  plots  of  ( 0 » 4* » 4s )  vs  6  ^or  <()=450  are  shown  in  Fig.  4-1  and 
Fig.  4-2.  The  spacing  for  which  each  detection  system  is  designed 
to  be  optimum  is  0.5a.  Also,  we  let  ea(t)=-eb(t)  so  that  the  binary 
signals  are  biphase  modulated.  Furthermore,  we  assume  that 

l  fT  ? 

(4-18)  l  e‘(t)dt  =  1 
£  M  3 

In  the  neighborhood  of  45°,  it  is  observed  that  |mLa(6*450;0) | 
increases  rapidly  and  eventually  reaches  the  maximum  at  45°. 
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These  two  figures  indicate  that  the  detection  system  is  capable  of 
suppressing  the  signal  a (t )  very  effectively  if  the  separation  of  o 
and  <j>  is  large.  As  the  incident  angles  e  and  o  get  very  close, 
the  suppression  becomes  very  poor. 

We  now  consider  the  case  when  the  signal  a(t)  propagates  into  the 
array  from  the  desired  angle  e.  The  ability  of  the  processor  to  favor 
this  signal  is  also  dependent  on  the  angular  separation.  In  Fig.  4-3 
and  Fig.  4-4,  we  plot  m|_a(0,<j>;e)  vs  e,  given  <p45°.  The  antenna 
spacing  used  by  each  system  is  again  equal  to  0.5x.  In  Fig.  4-3, 
where  m=4,  we  observe  that  |mj_a( e,450; e)  |  begins  to  decrease  steadily 
at  e=0°  and  0=70°.  In  Fig.  4-4,  where  m=12,  the  absolute  value  of 
each  processor  output  decreases  sharply  in  the  neighborhood  of  e=45° 
which  is  the  incident  angle  for  the  interfering  signal. 

For  a  constant  m  and  d,  the  difference  between  the  delay  times  t 
and  v  decreases  as  the  incident  angles  0  and  <j>  get  closer  and  closer. 
When  the  separation  between  0  and  </>  is  large  enough,  the  change  of 
the  difference  between  t  and  v  causes  the  sine  and  cosine  functions  in 
Eq.  (4-10)  and  Eq.  (4-15)  to  fluctuate  periodically.  When  0  and  <p  get 
closer  and  closer,  the  distinction  between  the  signals  from  e  and  <> 
diminishes  gradually.  Consequently,  the  processor  becomes  less 
capable  of  favoring  one  type  of  signal  and  suppressing  another  type 
and  vice  versa.  Finally,  the  two  signals  get  equal  treatment  when 
0=<fr. 


The  upper  bound  of  |mLa(0*45°;e) |  for  m=12  should  be  three  times 
higher  than  m=4.  This  property  is  reflected  in  Fig.  4-3  and  Fig.  4-4. 
One  of  the  differences  between  mLa(0»<J>;0)  and  mLa(e, <(>;<{>)  is  shown 
in  those  figures,  i.e.,  the  former  does  not  change  sign  whereas 
the  latter  does  change  sign  as  e  is  varied. 


We  pointed  out  earlier  that  m|_a ( 0, 4>)  vanishes  if  a, 4)  and  n 
satisfy  Eq.  (4-15).  It  can  be  shown  from  Eq.  (4-13)  that  Eq.  (4-15) 
is  also  the  condition  for  |nu_a( ©» 0)  |  to  reach  its  upper  bound. 

From  the  discussion  given  earlier,  it  appears  that  an  antenna  array 
with  spacing  d  satisfying  this  condition,  for  a  given  0,4,  and  m, 
would  optimize  the  error  rate  performance  of  the  system.  We  shall 
discuss  Eq.  (4-15)  more  later.  The  upper  bound  of  |m.  (e, 4>; <f>)  | 
increases  linearly  with  m.  The  function  ua 


(4-19) 


fOn.iy^1)  =  m  1 


which  appears  in  the  upper  bound  for  |mLa(0»<f>i<fr)  1 1  is  plotted  in 
Fig.  4-5  for  three  values  of  hw/hn.  It  is  seen  that  increasing  m 
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does  not  increase  ffm^h^  )  drastically.  We  also  observe  that 
fOn.hgh-1)  tends  to  approach  l^h"1  as  m  gets  large.  In  fact,  by 
expanding  Eq.  (4-19)  into  a  power  series,  we  find  that 

-1  hw 

(4-20)  lim  f(m,hwhn‘)  =  ^ 
m-**>  n 

Equation  (4-20)  suggests  that  the  upper  bound  for  I mi_a ( e , <t> ; )  |  can 
not  be  raised  appreciably  by  increasing  m.  Hence,  we  can  widen  the 
difference  between  mL  (e,(j>;e)  and  mL  (e,<r,4>)  by  using  a  large  array. 


20  24  28 


Fig.  4-5— f(m,hwh^  )  vs.  m. 


The  difference  in  the  magnitudes  of  |m,  (e,$;e)|  and  |m,  (e, <*>;$) 
can  be  best  indicated  by  the  ratio 


which  we  may  call  the  discrimination  figure  (D.F.)  because  it 
shows  the  extent  that  the  processor  favors  one  signal  and  suppresses 
another  signal.  The  D.F.  for  the  detector  using  a  4-element  array 
with  spacing  d  =  0.5x  is  shown  in  Fig.  4-6.  The  incident  angle  for 
the  Interfering  signal  is  45°.  As  can  be  seen,  the  D.F.  is  usually 
quite  high  if  the  separation  between  e  and  <f>  is  large,  but  decreases 
sharply  as  e  get  closer  to  4.  Ultimately,  this  figure  goes  to  zero 
at  0s4>,  because  m|_a(d,<t>;e)  and  m|_a(6,<t>;<j>)  become  equal  when  es<t>. 

For  a  constant  m,  the  D.F.  changes  with  the  antenna  spacing  d. 

Figure  4-7  shows  a  comparison  of  the  D.F.  for  d  =  0.5x  and  for 
d  ■  0.75X  given  m=2.  A  plot  of  this  figure,  given  e=60°  and  $=45°, 
of  a  4-element  array  detector  vs  tv/hn  is  shown  in  Fig.  4-8.  The 
quantities  |mi_a  ( 0 , 4> ;  e )  I  and  imi_a(0,<i);4>>  |  increase  with  hyjhn  at 
different  rates.  Since  the  denominator  in  the  expression  of  D.F. 
increases  faster  with  hy/ 1 hn  than  the  numerator,  the  D.F.  declines  as 
hwlhn  increases.  This  indicates  that  the  more  internal  noise,  the 
poorer  is  the  figure. 

To  broaden  the  study  of  the  directional  characteristics  of  the 
detection  system*  we  also  examine  the  behavior  of  m|_a(0,<|>;0)  at  an 
arbitrary  angle  e.  In  order  to  compare  easily  the  treatment  of  the 
signal „a(t)  from  0  with  that  of  the  same  signal  from  the  desired 
angle  e,  we  study  the  normalized  ratio  defined  as  follows: 

n»La(e,$;e) 

mLa(e,<t>;e) 

For  a  given  0  and  <p,  the  F(e,<j>;0)  vs.  0  plot  shows  relative  angular 
response  of  the  system  to  the  testing  signal  a(t)  from  an  arbitrary 
angle  0. 

Several  typical  curves  of  F(4O°,<f>;0)  vs.  0  of  a  4-element 
array  processor  are  shown  in  Figs.  4-9  through  4-12.  The  spacing 
d  for  the  array  is  0.5X  and  hw/hn  is  set  equal  to  0.2.  The  incident 
angle  <t>  for  thejnterference  in  Fig.  4-9  js  85°.  The  figure  shows 
that  F(4O°,8O°;0)  reaches  1$s  maximum  at  0=$  and  becomes  a  relative 
minimum  at  0=<|>.  But  F(e,$;e)  dges  not  always  behave  in  this  manner. 

In  Fig.  4-10,  where  <j>=45°,  the  0  which  equals  the  incident  angle  for 
the  desired  signal  does  not  correspond  to  the  maximum,  nor  does  that 
for  the  interference  correspond  to  a  relative  minimum.  Before  giving 
an  explanation  for  this,  we  show  the  response  curves  for  systems  using 
larger  antenna  arrays  in  Fig.  4-11  and  Fig.  4-12.  The  spacing  for 
each  system  Is  0.5x.  In  Fig.  4-11,  where  m=8,  the  discrimination 
figure  increases  tg  the  neighborhood  of  20  dB.  In  Fig.  4-12,  where 
m=12,  F(40°,  45°;  ©preaches  the  peak  value  near  0=0  and  becomes  a 
relative  minimum  at  0=$  even  though  the  separation  between  0  and  <p 
is  only  5°. 


F(0,<j>;0)  =  20*log 
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8  (DEGREES) 


Fig.  4-7~D1scr1m1  nation  figure  given  m=2,  <j>=45°  and  h  h'^0.2. 

Il  H 


for  all  e's.  From  Eq.  (4-12)  we  find  that  the  condition  given  above 
is  satisfied  if 


(4-21)  sin  5J-  w(t-v)  3  0 

where  t>*v(1.e.,e^).  It  is  interesting  to  point  out  that  this 
condition  ajso  makes  the  response  curve  become  the  absolute  minimum 
(-«  dB)  at  e=4>. 

Since  mia(e,*;e)  Is  equal  to  the  component  contributed  to  L(R) 
by  the  signal  a(t)  from  the  desired  source,  large  |m|_a (e ; e )  |  will 
reduce  the  influence  of  the  interference  on  making  decisions.  The 
fact  that  each  curve  in  Fig.  4-10  through  Fig.  4-12  fails  to  assume 
the  maximum  value  at  e*e  Indicates  that  the  condition  given  in 
Eq.  (21)  is  not  satisfied  for  each  corresponding  detector. 

The  optimum  system  Is  assumed  to  operate  with  one  desired  source 
and  one  interfering  source.  The  incident  angle  for  the  desired 
signal  is  e.  The  angular  response  curves  shown  in  Fig.  4-9  through 
Fig.  4-12  suggest  that  it  might  be  possible  for  a  certain  detector 
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to  be  optimum  for  two  desired  angles.  This  is  indeed  true  if  the 
second  desired  source  also  transmits  the  binary  signals  of  the  same 
waveforms.  This  will  be  shown  by  an  example  given  in  Section  IV. F. 

D.  The  Probability  of  Error 

The  discussion  so  far  concerns  the  nature  of  how  the  optimum 
detection  system  treats  the  signal  a(t)  from  different  incident 
angles.  With  this  background  we  next  examine  the  system  performance. 
For  a  digital  communication  system,  the  performance  of  the  system  is 
judged  mainly  by  its  ability  to  make  correct  decisions.  For  this 
reason,  the  technique  of  minimizing  Pe(np)  defined  in  Eq.  (3-32)  is 
investigated.  As  will  become  clear,  minimization  will  involve  choosing 
the  total  number  of  antennas  and  the  spacing  d  used  by  the  antenna 
array.  In  this  section,  we  assume  that  the  interfering  signal  is  a 
white  noise  process  with  spectral  height  h  .  The  desired  binary 
signals  are  assumed  biphase  modulated  in  tne  sense  that  ea(t)=-et,(t). 
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Fig.  4-ll~F(40°,45°;e)  vs.  0  given  m=8,  d=  ^  and  hyh^O.2 


Fig.  4-12— F(40o,45o;e)  vs.  0  given  m=12,  d=  j  and  h  h_1=0.2. 

^  W  fl 
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Hence,  the  difference  of  the  means  of  the  processor  output  L(R) 
under  hypotheses  Hb  and  Ha  is  equal  to 

f  1  sin2  J  o)(t-v)1  fT  Ceb(t)"ea(t)^ 


mLb‘‘n,La  " 


+  "w  sin*  £  (t-v) 


m  + 


(t-v)  fT 

-  . 

,-v)J  *T 


where 


j(t-v)  =  2iT  x  ^cos  0  _  C0S  ^ 


A  convenient  figure  for  showing  how  badly  the  interference 
and  noise  are  for  this  problem  is  the  ratio* 


which  in  fact  is  a  signal  energy  to  the  noise  power  density  ratio. 
By  making  use  of  the  assumption  given  in  Eq.  (4-19),  we  may  write 


This  equation  enables  us  to  specify  hy,  in  terms  of  the  spectral 
height  ratio  hw/hn  and  the  quantity  (S/N), 

We  now  consider  the  detection  system  which  usea  a  4-element 
array  and  is  designed  for  0=55°  and  <j>=45°.  The  array  spacing  d 
is  0.5x.  Table  4-1  shows  Pp,  Pq  and  Pe  vs  the  ratio  n/oi  (n  is 
the  threshold  setting).  This  system  is  assumed  to  operate  under 
the  conditions  that  (S/N)  =  1.78  (2.5  dB)  and  hn  is  five  times  higher 
than  1^.  The  spectral  height  for  the  internal  white  noise  is  assumed 
equal  to  0.0093.  From  Table  4-1  we  see  that  the  probability  of  making 
incorrect  decisions  assumes  its  lowest  value  0.6846  x  10"3  at »)  =0 
which  agrees  with  Eq.  (3-30)  because  mj_a  is  equal  to  -my}  here.  To 


The  notation  (S/N)  does  not  imply  the  ratio  is  the  signal-to- 
noise  ratio. 


V 


reduce  the  detection  error  P  (n0),  it  is  clear  from  Eq.  (3-31)  that 
we  should  increase  e 


The  value  for  (mj_b_mLg ) may  be  increased  by  adjusting  the  spacing 
d  or  by  varying  m  as  indicated  by  the  aquation  itself.  A  plot  of 
(mLb“mLa ) ab '  vs  m  1S  s^own  in  Pig*  4-13  for  three  different  values  of 
hw/nn.  The  angles  0  and  <j>  are  assumed  55°  and  45°,  respectively. 

In  plotting  this  figure,  we  also  let  (S/N)  =  1.78  and  d  =  0.5x. 

When  hw/hn  =  0.2,  (m^  -  mLa)a[‘  at  m=6  is  about  10.6,  which  is 
higher  than  the  value  at  m=4.  Hence,  we  expect  that  the  detection 
error  Pe(n0)  when  m=6  would  be  lower  than  0.6846  x  10~3.  in  fact, 
Pe(n0)  for  m=6  is  0.4640  x  10-7. 


TABLE  4-1 

PF,  PQ  and  Pr(e)  vs.  n/oL. 


n/oL 

PF 

PD 

p 

e 

-5.00 

.96388233E 

0 

. 1 0000000E 

1 

. 481 9E 

0 

-4.00 

.78745639E 

0 

. 1 0000000E 

1 

.3937E 

0 

-3.00 

.41981232E 

0 

. 1 0000000E 

1 

.2099E 

0 

-2.00 

.11460924E 

0 

. 1 0000000E 

1 

.5740E 

-1 

-1.00 

. 1 381 8979E 

-1 

.99998689E 

0 

. 691 6E 

-2 

0.00 

. 68460932E 

-3 

. 99931 550E 

0 

.6046E 

-3 

1.00 

.13207964E 

-4 

.98618102E 

0 

. 691 6E 

-2 

2.00 

. 9837951 OE 

-7 

.88539076E 

0 

.5730E 

-1 

3.00 

.27808902E 

-9 

.5801 8768E 

0 

.2099E 

0 

4.00 

.29586578E 

-12 

.21 254361 E 

0 

.3937E 

0 

5.00 

.1 1784304E 

-15 

.36117673E 

-1 

.4987E 

0 

The  internal  noise  degrades  the  error  rate  performance.  To 
show  this,  curves  of  Pe(on)  vs  hyy/h,,  for  the  detectors  (55°,  45°) 
with  m=4  and  6  respectively  are  shown  in  Fig.  4-14.  As  can  be  seen 
from  Eq.  (4-22),  the  spectral  height  hy*  increases  with  the  ratio 
hw/hn  if  the  quantity  (S/N)  is  constant.  In  plotting  Fig.  4-14  we 
have  let  (S/N)  =  1.78  which  is  about  2.5  dB.  In  Fig.  4-15,  we  show 
the  curves  of  Pe(n0)  vs  (S/N)  for  the  detector  (55°, 45°)  using  a 
4-element  array.  The  antenna  spacing  is  again  equal  to  0.5a.  It 
can  be  seen  from  these  curves  that  the  internal  noise  degrades  the 
error  rate  performance  drastically. 
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We  have  shown  that  it  is  necessary  to  increase  (mLb-mLa)°[ 
in  order  to  reduce  Pe(min).  For  the  example  given  above,  we 
accomplished  this  by  increasing  the  total  number  of  antennas.  But 
increasing  m  is  not  the  only  way.  As  a  matter  of  fact,  the  best 
approach  to  this  problem  is  to  choose  the  proper  spacing  for  the 
elements  of  the  antenna  array. 

E.  The  Optimum  Spacing 

It  has  been  pointed  out  earlier  that  for  the  case  when  the 
interference  is  a  Gaussian  white  noise  process  jm[_a  |  becomes  equal 
to  its  upper  bound  if  e,4>  ,  m  and  d  are  related  by  Eq.  (4-".5). 
Solving  Eq.  (4-15)  for  d,  one  finds  the  solution  is  many-valued 
because  k  can  be  any  integer.  We  define  the  optimum  spacing  dop 
as  the  solution  when  k=l  is  chosen,  i.e.. 


(4-24)  d._  A  7T— — — | ,  cos  e  -  cos  <f>  f  0. 

op  —  mjcos  e-  cos  <j>|  r 

Since  Eq.  (4-15)  is  not  valid  if  e=<j>,  there  does  not  exist  an  optimum 
spacing  when  e=«t>. 

Equation  (3-32)  indicates  that  the  error  Pe(no)  decreases  as 
(mLb“mLa)  *  fyy  *ncreases»  regardless  of  the  power  spectrum  of  the 
interfering  signal  or  the  waveforms  for  the  desired  signals.  When 
the  interference  is  a  white  noise  process,  we  conclude  from  Eq.  (4-23) 
that  the  smallest  value  for  Pe(oo)  "is  obtained  if  d  is  equal  to  the 
optimum  spacing.  The  reason  is  that  the  quantity  (m.  .-m,  Jar'  has  the 
peak  value  at  this  particular  spacing.  LD  La  L 

Referring  to  Eq.  (3-14),  we  see  that  the  amount  of  time-shifting 
on  s(t)  and  the  amount  of  time-delaying  on  r(t)  depend  on  the  incident 
angles  and  the  antenna  spacing.  When  the  spacing  is  changed,  we  must 
adjust  each  delay  line  to  provide  the  correct  delay. 

Referring  to  Eq.  (4-23),  we  find  that  (mLb-mLaW1  may  be 
rewritten  as 


where  the  ratio  d/dop  may  be  called  the  normalized  spacing.  In 
order  to  compare  the  error  rate  performance  for  systems  optimized 
for  different  spacings,  several  curves  of  Pe(n0)  vs.  d/dop  are 
shown  in  Fig.  4-16.  In  plotting  these  curves,  we  let  (S/N)  =  1.78 


- 1 


i 


Q 
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and  hw/hn  =  0.2.  In  Fig.  4-16,  each  curve  reaches  the  minimum  at 
d/d0n  =  1  as  we  have  expected.  The  error  Pe(no)  at  d/dop  =  1  for 
m  =  3  is  almost  one  thousand  times  lower  than  the  value  at  0.5. 

This  clearly  points  out  that  the  spacing  d  plays  a  very  important 
role  in  the  error  rate  performance.  We  shall  use  Pe(min)  to  denote 
the  error  rate  for  the  system  which  is  designed  for  an  optimum  spacing. 

To  show  the  lowest  detection  error  available  to  a  given  array 
size,  the  curves  of  P£(min)  vs.  (S/N)  for  m=2  through  5  given 
h^/hn  =  0.2,  0.5  and  1.0  are  shown  in  Fig.  4-17  through  Fig.  4-19. 

These  curves  are  valid  for  arbitrary  incident  angles  0  and  <j>  as  long 
as  the  spacing  used  by  each  system  is  equal  to  the  optimum  spacing 
for  the  given  e  and  <j>.  The  reason  for  this  is  given  in  the  next 
section. 


F.  Arbitrary  Power  Spectrum 

The  optimum  spacing  d0p  defingd  by  Eq.  (4-24)  is  found  by 
examining  the  behavior  of  mLa(0»<Pi0)  under  the  assumption  that  the 
interfering  signal  is  a  white  noise  process.  In  this  section,  we 
shall  show  that  this  particular  spacing  also  optimizes  the  error 
rate  performance  even  though  the  interfering  signal  is  characterized 
by  a  different  power  spectrum.  The  signal  waveforms  considered  in 
this  section  are  assumed  to  have  the  general  expressions  as  given 
in  Eq.  (4-1)  and  Eq.  (4-2). 


For  an  arbitrary  power  spectrum,  the  difference  between  the  means 
of  the  test  function  L(R)  under  hypotheses  %  and  Ha  may  be  derived 
from  Eq.  (3-14).  The  result  may  be  written  as 

(4-25) 


Lb  •  mLa  AE(L(R)|Hb}  -  E TL(R)  | Ha } 
m  rT  9  m  f  rT 

=  l  s?(t-)dt  -  [n  s.(t.)s,(z.)*')'(t,z)-dt  dz 

i=l  j-T  1  1  i,j=lJj-T  1  1  J  J 


=  I 


m  rT  9  00  1  fT  m 

I  s?(t.)dt  -  l  — \ — S  l  si(t.)<t»  (t)dt 

1*1  J-T  11  <1=1  m  .  hw  J-T  i=l  1  1  4 

m  +  t — c— 


w 


Since  the  eigenvalues  ( Ag,-hw)  are  always  positive  and  s(t)  and  <f>^(t ) 
are  real  functions,  we  may  write  that 


(4-26) 


m 


Lb  '  mLa 


m 


-T 


r(ti )dt 


53 


f 


1 


The  closed  form  expression  for 


I  sf(t.) 
i=l  1  1 

is  readily  obtained  when  we  replace  t  and  ea(t)  in  Eq.  (4-10)  by  T 
and  e(t),  respectively.  Hence, 


(4-27) 


m  rT  rT  2 1.\ 

l  s?(t.)dt  =  m  •  f 

i=l  J-T  11  J-T  6 


» 

in  maj( x-v)  fT  e2(t)  /I  m-1  \ 

Tinsn^r  *  j_T  ~z~  * co^2^  -  -r 

+  2<(>s  -  (m-1  )urrj’  dt 


The  second  term  of  Eq.  (4-27)  is  much  smaller  than  the  first  term. 
Thus  the  upper  bound  of  (mLb-mLa)  is  approximately  equal  to 


fT  2 

e£(t)  •  dt. 
J-T 


When  the  array  is  spaced  according  to  Eq.  (4-24),  we  find  that 


(4-28)  sin  Sy-  o>(t-v)  =  sin  mw(T-v)  =  0 


As  a  result,  we  conclude  from  Eq.  (4-8)  that 


i  MV  *  o 

i=l  1  1 

Consequently,  the  left  side  of  Eq.  (4-26)  is  equal  to  the  right  side. 
Also,  the  second  term  of  Eq.  (4-27)  vanishes  because  sin  mu)(-r-v)  =  0. 
Hence,  we  obtain 


m  rT  9 

(4-29)  mLb  -  \a  =  7  eZ(t)  dt 
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which  is  very  close  to  the  highest  value  of  (mLt>-m|_a)  we  can  get. 
Substituting  Eq.  (4-29)  into  Eq.  (3-32),  we  obtain  che  detection 
error 


(4-30)  P£(min)  =  erfc*^Jjj-  j  ^  dt^ 


This  is  the  lowest  error  rate  available  for  a  given  array  size. 

Equation  (4-30)  does  not  contain  the  incident  angles  0  and  4>. 
This  implies  that  this  error  rate  does  not  depend  on  the  angular 
separation  between  the  interfering  source  and  the  desired  signal 
source.  The  optimum  spacing  is  a  function  of  e,  <j>  and  m.  To 
obtain  this  lowest  error  rate,  we  need  to  change  the  antenna  spacing 
and  adjust  the  delay  lines  when  the  incident  angles  e  and  <p  vary. 

We  also  notice  that  Eq.  (4-30)  does  not  contain  the  parameter 
which  characterizes  the  interfering  signal.  This  points  out  that 
this  lowest  possible  error  rate  is  independent  of  the  power  spectrum 
of  the  interfering  signal.  The  reason  is  that  at  this  particular 
antenna  spacing  the  interfering  signal  is  completely  cancelled  by  the 
optimum  processor.  A  formal  proof  of  this  statement  will  be  given 
in  Chapter  VI. 

In  plotting  the  curves  in  Fig.  4-17  through  Fig.  4-19,  we  have 
used  the  optimum  spacing  for  the  antenna  elements.  Hence,  these 
curves  are  also  valid  for  other  types  of  power  spectrums  if  the 
desired  signals  are  biphase  modulated  in  the  sense  that  e(t)=-e.(t). 

a  D 

If  we  let  the  signal  a(t)  =  0,  then  Eq.  (4-1)  and  Eq.  (4-2) 
represent  an  on-off  keying  binary  signals.  The  structure  for  the 
optimum  processor  becomes  simpler  as  can  be  seen  by  the  test 
statistics 

(4-31) 

m 

L0(R)  A  I 

0  1*1 


•  v(t,z)  •  dt  dz, 

which  is  obtained  by  replacing  s(t)  in  Eq.  (3-14)  by  b(t).  Ac¬ 
cordingly,  the  threshold  that  minimizes  the  detection  error  is 

\  *  7  mLb 


[T  -  n  j|T  • 

J_T  '»j-1  _T 
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The  lowest  error  rate  for  a  given  array  size  may  be  derived  from 
Eq.  (4-30).  We  find  that 


(4-32) 


P£(min)  = 


erfc. 


f1  *  fT  e>(t)  dt') 

UK  J.t  -T-  dt; 


Hence,  for  this  on-off  binary  signals  detector,  the  error  rate 
performance  can  be  made  as  good  as  the  biphase-binary  detector  by 
increasing  the  amplitude  of  the  signal  b(t). 

It  was  pointed  out  in  Section  IVC  that  some  detectors  may  be 
optimum  for  two  desired  sources  if  they  transmit  the  same  type  of 
binary  signals.  We  now  consider  a  detection  system  which  is 
designed  to  be  optimum  for  e=cos“'  (0.457),  <p=45°,  m=4  and  d/A  =  1. 
We  find  that  the  conditions  given  in  Eq.  (4-28)  are  satisfied  by 
the  data  given  above.  Hence  we  obtain 


(4-33) 


m, 


Lb 


-  m 


La 


=  2 


-T 


e^(t)dt 


When  the  same  binary  signals  propagate  into  the  array  for  the  same 
detection  systems  from  angle  e=cos-l  (0.957),  we  find  that  Eq.  (4-28) 
is  again  satisfied.  Hence,  the  difference  of  the  means  is  identical 
to  Eq.  (4-33).  Apparently,  the  same  system  provides  the  same  error 
rate  performance  to  two  incident  angles. 
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CHAPTER  V 

THE  DETECTOR  FOR  BINARY  SIGNALS  HAVING 
A  RANDOM  PHASE  AND  AMPLITUDE 


A.  The  Optimum  Processor 

This  chapter  considers  the  case  where  the  desired  signal  is  no 
longer  assumed  to  be  known  exactly,  but  instead  has  certain  random 
parameters. 


In  a  realistic  situation,  there  are  always  at  least  some 
parameters  characterizing  the  communication  signal  that  are  not  known. 
For  instance,  in  high  frequency  communi cation  via  the  ionosphere  and 
in  channels  employing  tropospheric  propagation,  the  signal  received 
usually  fluctuates  in  amplitude  and  phase  due  to  fading  and  multipath 
effects.  These  effects  are  usually  modeled  by  assuming  certain 
parameters  of  the  communication  signal  to  be  random  variables.  Signal 
detection  in  this  case  is  often  referred  to  as  the  composite  hypothesis 
-  testing  problem. 


The  technique  for  determining  an  optimum  processor  for  detecting 
the  binary  signals  having  random  parameters  is  straightforward,  because 
the  only  new  feature  is  finding  p( R | H^) ,  where  k=a  or  b,  when  certain 
parameters  are  random.  For  the  case  when  the  random  parameters  have 
the  same  statistics  on  both  hypotheses,  the  likelihood  ratio  is  simply 
given  by 


(5-1) 

L.R. 


P(R|Hb) 

P(ft|Ha) 


p(R|r,H.)  •  p(r)  •  dr 

=  lim  - ° - 

N-*»  p(Rjr,H  )  •  p(r)  •  dr 

C 


where  p(r)  is  the  probability  density  function  for  the  parameter 
vector  r,  and  p(R|r,l!k)  is  the  density  function  of  R  under 
hypothesis  Hk  for  a  given  r  [22], 

In  this  chapter  we  consider  the  binary  signals  of  the  form: 

Hfl:  a(t)  -  f  •  efl(t)  •  cos(u>at  +  $s) 

(5-2) 

Hb:  b(t)  =  f  •  eb(t)  •  cos(wbt  +  <f>s) 
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where  the  amplitude  f  is  a  Rayleigh  random  variable  and  the  phase 
<j>$  is  assumed  uniformly  distributed  between  -ir  and  it. 

Three  comnonly  used  methods  for  transmitting  binary  signals  over 
an  additive  Gaussian  noise  channel  are  on-off  keying  (ASK),  frequency- 
shift  keying  (FSK)  and  phase-shift  keying  (PSK).  It  is  apparent  that 
the  model  we  have  here  can  be  applied  to  either  the  ASK  or  FSK  case. 

The  signals  given  in  Eq.  (5-2)  may  be  rewritten  as 


a  ( t ) 

r 

03 
73 
i _ 

va(t)  ' 

=  x  • 

-  y  • 

.b(t). 

-  vb(t>- 

where  we  define 


(5-4) 

X  4  f  •  COS  <j>s 

(5-5) 

y  4  f  •  sin  4>s 

(5-6) 

uk(t)  ek(t)  •  cos  oikt,  k  =  a 

and  b. 

and 

(5-7) 

Vk(t)  a  ek(t)  •  sin  wkt  . 

x  and 
equal 

y  are  independent  zero-mean  Gaussian 
variances,  say  aj!  C23].  Let 

random  variables  with 

then  the  density  function  for  the  random  vector  r  may  be  written 
as 


(5-8) 


P(r)  = 


Since  u|<(t)  and  V|<(t)  do  not  contain  unknown  parameters,  they  are 
completely  known  functions.  Moreover,  they  are  90°  out  of  phase. 
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The  desired  signal  may  then  be  considered  as  the  sum  of  two  signals 
for  which  the  amplitudes  are  Gaussian  random  variables. 

Under  hypothesis  the  density  of  R  associated  with  an 
m-element  antenna  array  and  an  N-term  Karhunen-Loeve  expansion 
is  equal  to  ?  ? 

-^r 

^  ff  2ar 

(5-9)  p(R|H.  )  =  -U-  p(R|r,Hh)  •  e  •  dx  dy 

D  2iro*  JJ  D 

r  -°° 

In  Eq.  (5-9),  the  conditional  density  is  given  by 


(5-10) 


expl-  1  (R-B)T  .  Q  .  (R-B)} 
|r,H.  )  =  lim  —h==— - *- 


N"~  J(2n)mN  | K | 


=  lim  — - 

l(27r)mN|K| 


f  . 

exp-j-  i  I  I  ff  rJt)rAz)  •  H?.(t,z)  •  dt  dz 
Z  i,j=l  jj-T  1  J  1J 


m  rT 

+  II  r.(t)b.(z)  •  H?.(t,z)  •  dt  dz 

i,j=l  1  J  1J 

,  m  r  rT  'j 

-  v  II  b. (t)b.(z)  •  H?.(t,z)  •  dt  dzr 
i  ,0=1  Jj-T  1  3  1J  J 


From  Eq.  (5-3),  we  obtain 
(5-11) 

t>,(t)b.|(z)  =  x2  .  ub1(t)  ubj-(z)  *  y2  •  vb1  (t)  vbj-(z) 

-  xy.[ub1(t)  ubJ(z)  +  vo1(t)  »bj(z)] 

In  terms  of  the  in-phase  and  quadrature  components  of  b(t),  we  may 
rewri te 
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where  we  define 


(5-13) 


m 


i,j=l  "-T 


%  A°r  t  ubi(t)  •  ub.(z)  •  H?j (t,z)  •  dt  dz 


bj 


1J 


2 

%  1 


(5-14) 


i.J 


(5-15) 


m  rT 

p 

m  rrl 

I 

1-1  J- 

r  ubA->dt- 

II 

i,j=l  JJ. 

I  •  dt 

4 

,  III 

r  ybi(t>  •  vm 

U)  .  H"j( 

f. 

rT  ? 

m 

U  • 

LTVbi(ti)dt  " 
\ 

II 

i »j=l 

.z)  • 

dt  dzi 

m 


Yb  -  °r  ii 


i ,  J=1 


ubi(t)vbj(z)  *  H^(t,z)  •  dt  dz 
Op  m  rT  m  f  rT 

’  ^  J-T  Ubi(ti)vbi(t1,dt  '  ii,  }{.><», >VZJ> 

•  v( t , z )  •  dt  dzl 


Lbu(r)  &  H  JJ  •  H.  ^ (t,z)  •  dt  dz 

i »j-l  ~T 

(j,  fT  ri(ti)ubi<ti)dt 

m  rfl 

I  I  r . (1; - ) *u.  . (z - )  •  *(t#z)  .  dt  dz  ' 

i,j=l  JJ-T  1  1  DJ  3 


(5-17) 


Lbv(r)  A  -  II  |J  rl(t)vbj(z)  *  *  dt  dz 

°~k  '  1=1  l!T  r1(t1)VM<t1)dt 

m  rrT 

-  II  *  'i'(t»z)  ‘  dt  dz  ' 

i ,j=l  JJ-T  1  1  DJ  J 


(5-18)  L  (r)  A  -  I  H  f|  r - (t )  r.(z)  •  H?.(t,z)  dt  dz 
r  6  i,j=l  JJ-T  1  3  13 


Substituting  Eq.  (5-12)  into  Eq.  (5-9),  yields 


(5-19) 


Lr(r) 


p(R|H.)  =  — U-  lim  e. -  = 

27ror  J( 2Tr)mN  |  k| 


exp  1-  j-?  x2(ab+l)  +  yt(6b+l)-2xy  •  Yb 


-  2xo2  L  Jr)  -  2y  a2 


k„(r)"|i  • 


dx  dy 


The  parameters  ab»  6b  anc*  Yb  are  not  random  variables.  For  a 
given  interfering  power  spectrum  and  variance  a2,  their  values  are 
controlled  by  the  spacing  d,  spectral  height  hw,  number  of  element 
m,  etc.  Lbu(r)  and  LkV(r)  are  Gaussian  random  variables,  containing 
the  information  needed  for  decision  making.  From  their  expressions, 
we  see  that  Lbu(r)  and  Lby(r)  may  be  realized  by  schemes  similar  to 
the  one  shown  in  Fig.  3-1.  It  will  become  clear  later  that  the 
random  variable  Lr(r)  is  irrelevant  in  the  hypothesis  testing. 
Equation  (5-19)  may  be  simplified  by  employing 


,oo 

. 

—  00 


exp< 


dz  =  2 rra 


The  procedure  for  simplification  is  straightforward  but  tedious.  Let 


(5-20)  Ab(Lu,Lv)  a  (a5+^)*Lbv(r)  +  2Yb'Lbu(r)Lbv(r J  +  b+1  >  *  Lbu ( r) 

and 

(5-21 )  u  b(a  ,6  »y  )  A  b+^  *  (&b  +  1 )  “  Y  b  * 


We  find  that 


(5-22) 


L>) 


p(R|  H.)  =  lim 

D  it 


ur 

T 


Ab^Lu*Lv^ 

P  b(a  ,6  <Y  ) 


>  (2n  )mN|  K  |  JT^TT^Y 


Replacing  the  subscript  b  by  a  in  Eq.  (5-13)  through  Eq.  (5-17),  we 
generate  a  set  of  new  functions  which  are  the  counterparts  of  aj,,  $b> 
Yb»  l-bu(r)  and  Lbv(r)  under  hypothesis  Ha.  We  denote  them  byaa,  Ba, 
y a*  Lau (r )  ancl  Lav(r)»  respectively.  For  example 


m  ffT 

jljZ,  Jj_T  ua1(t)  •  uaj(z>  ‘  *  dt  dz 


The  density  of  R  under  Ha  may  be  formulated  in  the  same  manner. 
p(R|Ha)  is  readily  obtained  if  Ab(Lu.Lv)  andpb(a,  6*y)  in  Eq.  (5-22) 
are  replaced  by  their  counterparts  a a(Lu»Lv)  and  v a(a,3,Y),  etc. 
Hence,  the  likelihood  ratio  is  given  by 


(5-23) 


L.R.  =  lim 


p(R|hk) 


N-x»  p  (R  |H_) 

a 

Ma(a,3,y) 


f°?  ak(l,..lu) 


Ma(a,3,y)  ^or  Ab(Lu,Lv)  Aa(Lu>Lv) 

Pb(a»3»y)  exP  ’  Y~  "  vib(a,8,y)  ?  *  ua(a,6,y) 


Taking  the  logarithm  and  incorporating  common  terms  and  known 
functions  in  the  threshold,  we  obtain  the  equivalent  test 


(5-24) 


A q(r)  =  na(<*,B,y)  •  Ab(Lu,Lv)  -  yb(a,8,y)  •  Aa0-U.LV) 


<  threshold. 

Ha 

The  expression  for  the  test  function  Ag(r)  is  quite  complicated. 

The  decision  process  relies  on  observing  four  sufficient  statistics 
Lku(r)  and  Lkv(r)»  k  =  a  and  b,  which  are  Gaussian  random  variables. 
From  Eq.  (5-24),  we  note  that  the  optimum  processor,  in  general, 
consists  of  two  similar  portions,  i.e.,  Aa(Lu,Ly)  and  Ab(Lu,Lv). 

A  scheme  for  generating  AbUu»l-v)  is  showri  in  Fig.  5-1. 

B.  The  Distribution  function 

To  find  the  statistics  of  the  test  function  Ag(r),  we  need  to 
determine  its  distribution.  In  theory,  the  distribution  of  A  (r) 
under  each  hypothesis  may  be  calculated  by  the  integral  9 

(5-25) 

FA(AlHk)  =  p(Lbu,Lbv»Lau,LavlHk)  *  dLbu*dLbv’dLau*dLav 


where  ^is  the  hyperplane  defined  by  Eq.  (5-24)  and  p(Lbu»Lbv»Lau»Lav) 
is  the  joint  density  of  the  random  variables  L...(r),  L_..(r),  L.  7r) 
and  Lby(r)  under  hypothesis  Hk.  au  av  Du 


The  statistics  of  L|cu(r)  and  L|<v(r),  where  k=a  and  b,  under 
both  hypothesis  are  given  in  Appendix  B.  In  general,  they  are  cor¬ 
related  Gaussian  random  variables  with  zero-means.  Thus  the  4  x  4 
covariance  matrix  associated  with  the  density  function  p  (Lb  ,  Lb  , 
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(  ab  +  1 ) 


(  3b+  1 » 


Fig.  5-1 --A  scheme  for  generating  Ab(Lu»Lv). 

Lau*  Lav)  does  not  have  an>  zero  element.  We  expect  that,  in  general, 
the  inverse  of  this  covariance  matrix  will  not  be  diagonal.  Con¬ 
sequently,  the  expression  for  the  density  function  is  quite  com¬ 
plicated.  As  can  be  seen  from  Eq.  (5-24),  the  geometry  of  the 
hyperplane  is  not  simple  either.  For  these  reasons,  we  can  see  that 
the  quadruple  integral  in  Eq.  (5-25)  is  difficult  to  evaluate. 

For  a  constant  Ag(r),  Eq.  (5-24)  is  the  equation  of  an  ellipsoid 
of  four  dimensions.  To  simplify  the  geometry  ofT*,  we  shall 
diagonalize  the  quadratic  equation  defining n,  The  test  function 
a  (r)  may  be  written  as 


Vr)  = 


T  r 

Lbv(r>' 

Lbu(r) 

Lav(r) 

Vr). 

vH+1)  vYb 


rYb 

0 

0 


0  0 

ua(eb+U  o  o 

0  v(aa+1)  VYa 

0  VYa  V(8a 


Lbv(r> 

Lbu(r> 

Lav'r> 

Lau<r> 
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Apparently,  the  coordinates  which  diagonalize 


(5-26) 


Vr)  A 


- 1 

r~ 

7T 

< 

-s 

_ 1 

T  ' 

[_  Lku(r)  . 

VK+1)  VYk 
VYk  pq‘^6k+1 


Lkv(r> 

Lku(r>. 

where 


and 


k  ^  q 

k  =  a  and  b 

q  =  a  and  b. 


will  also  diagonalize  the  quadratic  equation  Ag(r). 


Let* 


(5-27) 


°k  A] 


o  u 


A  ^ 


k2-J 


“k  +  6k  +  2  *  6k 


ak  +  \  +  2  +  6k 


,  k  t  q. 


where 


«L  A  J(a.  -  8.  )2  +  4y l 


Jk  -  ^'  k  Pk' 

Also,  let  us  define 


(5-28)  G.  4  - 

/5T 


J-vv*k 

A_ 

9kll 

9kl2 " 

Lj»k-vsk 

J"“k+6k*sk  J 

-9k21 

9k22 . 

* 

Ckl  and  Ck2  are  the  eigenvalues  associated  with  the  2  x  2  symmetric 
matrix  given  in  the  right  side  of  Eq.  (5-26). 
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G|<  is  an  orthonormal  matrix.  Then  we  find  that  in  terms  of  the 
normal  coordinates  ofbu(r),cCbv(r),o^au(r)  andjfav(r),  where 


4v<r> 

"Lkv(r)' 

(5-29) 

=  Dk  *  Gk  * 

■  4,<r>- 

-Lku(r>- 

the  test  function  Ag(r)  has  the  simple  representation, 

Hb 

(5-30)  Ag(r)  =4v(r)  ^(r)  -£\v{r)  -^|u(r)  5  threshold. 

d 

As  a  result,  the  distribution  of  A  (r)  may  be  evaluated  alternatively 
by  9 


r  r  r  r 


R 


The  statistics  of  L|<u  and  L^y,  k=a  and  b,  under  both  hypotheses 
are  given  in  Appendix  B.  In  general,  they  are  correlated  Gaussian 
random  variables  with  zero-means.  Being  linear  functions  of  Lbu» 
Lbv»  Lau  and  Lav,  the  new  coordinates  are  therefore  also  zero-mean 
Gaussian  random  variables.  The  covariance  matrix  associated 
with  the  normal  coordinates  is  a  4  x  4  matrix  * 


(5-31) 


V 

^bu 

^buLbv 

^buAiu 

> 

o^bv*"bu 

i2 

*-bv 

^bv4u 

^bv^av 

^auLbu 

^au  bv 

^2 

dau 

^ai/av 

,'l-avLbu 

^av  bv 

4v4u 

^av 

From  Eq.  (5-29)  we  find  that  in  general  only 


(5-32)  B^u^vl  =  °.  k  *  a,  b. 
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For  this  reason  the  four  jointly  normal  random  variables i/s  cannot 
be  statistically  independent,  The  variables,^ 's  in  the  expression 
for  the  joint  density  p(i!bu»^bv»^au»iav)  are  only  partially  de¬ 
coupled.  As  a  result  of  this,  the  quadruple  integral  Eq.  (5-31) 
remains  difficult  to  integrate.  It  appears  that  a  simple,  analytical 
formula  for  the  distribution  of  the  test  function  A  (r)  given  by 
Eq.  (5-24)  seems  not  obtainable.  9 


C.  The  Probability  Equations,  ASK 

One  commonly  used  method  for  transmitting  the  binary  symbols  0 
and  1  using  ASK  signals  is  to  send  a  specific  signal  for  the  symbol 
1  and  to  send  no  signal  at  all  for  the  symbol  0.  For  this  special 
case,  only  two  test  variables  remain  in  the  equivalent  likelihood 
ratio  test  defined  by  Eq.  (5-24).  Besides,  the  distribution  of  the 
test  function  may  be  obtained  by  evaluating  a  double  integral. 

Suppose  the  signal  b(t)  is  sent  for  the  symbol  1.  The  signals 
received  by  the  ith  element  under  both  hypotheses  are 

H  :  r.(t)  =  n.(t)  +  w.(t) 

(5-33)  a  1  1  1 

Hfa:  ri (t)  =  ni(t)  +  w.(t)  +  x*ubi(t)  -  y  •  vbi(t) 

Setting  u,(t)  =  v_(t)  =  0,  we  obtain  from  Eq.  (5-24)  that 

a  a 

(5-34) 

A(r)  4  (y-1)  •  lby(r)  +  2i,b  •  Lbv(r)Lbu(r)  +  (Bb+1)  •  L*u(r) 


The  distribution  of  A(r)  may  be  evaluated  by 

FA(»|Hk)  '  }}  P«U^ulHk)  ‘  '  «4U 

We  shall  see  later  that  even  for  this  special  case,  the  double 
integral  given  above  can  not  be  evaluated  analytically  either. 
However,  comparatively  speaking,  the  expression  for  the  distribution 
is  quite  simple. 


The  structure  of  the  optimum  processor  for  detecting  on-off 
binary  signals  is  just  half  the  general  system  In  the  l_bu(r) 

and  Lbv( ^  plane  the  hypothesis  Ha  is  chosen  whenever  A(r)  appears 
within  the  ellipse  or  Hb  is  chosen  whenever  otherwise.  The  geometry 
(e.g.,  major  axis)  of  the  ellipse  is  decided  by  the  threshold, a  b»  Bb 
and  Yb.  From  Eq.  (5-29)  we  know  that  the  quadratic  equation  Eq.  (5-34) 
may  be  rewritten  as  u 


(5-35)  A(r)  =«f^(r) 


where* 


(5-36)  Jfbv  j  •  (9bn  Lbv  +  9b12  Lbu^ 


(5-37)  «4u  -j?b2  •  (gb2]  Lfav  +  gb22  Lbu) 

Making  use  of  Eqs.  (5-36)  and  (5-37),  we  find  that  the  variances  of 
4v(r)  and*bu^r)  under  hypothesis  Hk,  k=a  and  b,  are  equal  to 

(5-38) 


i,kAE2(^y(r)|Hk) 

=  (“b +  h  +  V  *  [7  (2+vv6b)J 


(5-39) 


°buki  var  {4u(r)|Hk} 


’b2 


K  +  -  O 


T7  *  '“b  T  pb  "  V 

2o„ 


e 

[7  (2+VV6b)] 


♦Because  the  factor  ua(“»0»Y)  does  not  appear  in  Eq.  (5-^4),  the 
eigenvalues  Cbl  and  ^2  EP*  (5-36)  are  equal  to  H(ab+0b+2-<sb) 
and  »s(ab+8b+2+6b),  respectively. 
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where 


ek  *  0  if  k  =  a 


ek  =  1  if  k  =  b 


A(r)  is  nonnegative  because X by  and^fbu  are  real  functions.  The 
distribution  function  F A( A jH. )  is  simply  given  by 

K  2  2 
(5-40)  _ibv  _  ^  bu 

F,(A|HJ  =  - ff  e  2obvk  2obuk  .  <g[  di 

A  k  2™bvk0bukjJR  bv  ^bu 


2  2 
(z  sin  e)  (z  cos  e) 


*A  t  2n  2 of 


‘^bvk^uk  JQJe=0 


•  z  dz  •  de 


The  probability  density  function  of  the  test  function  A(r)  under 
hypothesis  thus  is 

(5-41) 


.  . . 3FA<AlHk>  _  1  3FaWM 

Pa^IV  —  3A - =  7ff“ - 


tirobvk  °buk  JO 


z=JX 


A | sin2e  +  cos2x 


°bvk  °buk 


where  a  ^  0.  This  integral  can  not  be  integrated  analytically.  One 
may  use  a  numerical  integration  technique  to  get  the  value.  Denoting 
the  probability  of  choosing  hypothesis  Hi.  when  Hb  is  present  by  Pd, 
and  the  probability  of  choosing  Hb  when  Ha  is  present  by  Pf  and 
referring  to  Fig.  5-2,  we  find  that  ' 


(5-42) 


n 


whose  expression  is  similar  to  Eq.  (5-42)  and  can  be  readily  obtained 
from  the  same  equation  when  the  third  subscripts  of  abub  and  abvb 
are  changed  from  b  to  a.  Before  presenting  some  numerical  results, 
we  shall  briefly  discuss  the  parameters  ak,  Bk  and 

D.  On  Parameters  a^,  ek  and  yk 

Examining  the  equations  for  Pd  and  Pf,  we  note  that  they  are 
mainly  controlled  by  the  parameters  otb,  Bb  and  yb.  In  the  case 
of  the  general  optimum  detection  system  Eq.  (5-24),  we  expect  that 
the  problem  of  optimizing  the  system  performance  would  become  a 
problem  of  designing  proper  values  for  ak,  Bk  and  yk,  where  k=a  and 
b.  For  this  reason,  we  study  the  properties  of  parameters  a^»  Bb 
and  yb.  In  practice,  the  carrier  frequency  u>b  is  much  greater  than 
the  one  used  in  the  modulating  signal  eb(t).  Very  little  error 
would  be  introduced  when  we  approximate  ubl- (t.j)  and  vbl- (t^ )  by 

(5-44)  ubi (ti ^  *  eb^  *  cos{(ob[t-(m-l)v>(i-l)u>b(T-v)} 
and 


(5-45)  vb- (ti )  £  eb(t)  •  sin{o)[t-(m-l)v>(i-l)ajb(T-v)} 
Thus  we  obtain  the  following  closed  forms: 
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(5-46) 


m 


J,  ubi(ti) 


m 

s  i  n  ■x-  a)L 


7  “b 


(t-v) 


sin  -j—  (t-v) 


*  6h(t)*cos{wb!It  -  (t+v)3) 


(5-47) 


Jl  Vb1(t1> 


c'n  7  U)b^T"v) 
wb 

sin  (x-v) 


6[J(t)*sin{w|3Ct  -  (t+v)]} 


and 


(5-48) 


m 


sin  mu^Cx-v)  eb(t) 


i|1  ubi(Vvbi(V  =  s'in  u^fx-v)  *  “2  sin{2u,b[t  -  ^  (x+v)]} 


Comparing  their  definitions  with  the  mean  of  L(R),  Eq.  (3-4),  we  con¬ 
clude  that  parameters  ab,  6b  and  mLa(e,<j>;e)  possess  the  same  di¬ 
rectional  property.  Using  the  results  just  obtained  for 


m  m 

ubi^V*  i|1  vbi  (ti  ^ 


and 


m 

ubi^V  vbi^i^ 

we  may  derive  the  closed  form  representations  of  ctb,  3b  and  Yb* 
For  the  sake  of  saving  space  we  write  them  in  a  corrmon  expression 
designated  as  g(e-j »e25- 
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(5-49) 


g(e1»02)  A  *  111  *  sinteg-e^  • 


T  ent) 


-T 


T 


dt 


2  2 
oj;  sin  majj3(T-v)  rT  eb(t) 

+  *  sin  a>b{x-v)  •  j_T  -T-  • 


sin{2w^[t  -  (t+v)]  +  Oi+e^  1  •  dt 


sin  £  wb(-r-v) 

w. 

sin  y~  (t-v) 


r 

V 


cos 


-T 


eb(t)ebU)  * 


ft  -  (t+v)  + 


m-1 


•  sin{a)b  [z  — (t+v)]  +  e2}  •  v(t,z)  •  dt  dz 

where  e-|  and  e2  assume  different  values  under  different  circumstances. 
In  terms  of  g(e-|,e2),  we  find  that 


(5-50)  «b  =  g(0  ,  tt/2) 


and 


Bb  =  g(-  tt/2,  0) 


Yb  =  g(o,o). 


The  function  ub(t)  is  real,  so 
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m  f  (i 

i  J J_t  Ubl(V  ‘  ubj(zJ>  ' 


»  i  l  I  uk,(t,.)<j>„(t)  •  dt 

-  I  itL 

*=1 

m 


From  Eqs.  (5-13)  and  (5-14),  we  thus  establish  that 

(5-5,)  |  f  •£,<*,>« 

D  nw  i=l  j-T 


or  m  [1  0 

(5-52) 


The  upper  bounds  in  Eq.  (5-51)  and  Eq.  (5-52)  are  nonnegative  because 
the  integrands  ut^ti)  and  vg^tj)  are  nonnegative.  As  a  matter  of 
fact,  these  bounds  are  just  equal  to  the  sum  of  the  first  two  terms 
of  Eq.  (5-49).  For  the  case  when  the  interference  is  a  white  noise 
process  with  flat  spectrum  height  hn,  Eq.  (5-49)  reduces  to 


(5-53 ) 


V  fT  ef(t) 

gw(0i»e2)  =  TT  m  *  s1n(e2"6^  *  -g-  d 


Op  sin  mw^T-v) 
\  *  "sin  u>b(T-v) 


T  e‘(t) 


•  sin{2a)^[t  -  (t+v)3+0i  +  ©2 )  dt 

4  i  f  sin  %.(,-»)  I2 


+  HtI  s 


(5-53) 

(Cont.) 


Discarding  numerically  small  terms  in  Eq.  (5-53),  we  obtain  the 
following  approximations  good  for  the  white  noise  interference: 


When  the  spacing  d  between  two  adjacent  antennas  is 


where  is  the  wavelength  corresponding  to  the  angular  frequency  u>h, 
then  D 


(5-56)  sin  Sy-  ^(t-v)  =  sin  mos^T-v)  =  0 

The  spacing  d0p  has  been  referred  to  as  the  optimum  spacing  for  the 
given  e,  <f>  and  m.  As  a  result  of  placing  an  array  at  this  particular 
spacing,  we  find  that* 


m  f\,  m  ^ 

(5-57)  “bi^V  =  i|1  vbi^V  =  0 


•ft 

The  definition  of  ti  is:  tj  a_  t-(m-i)v  =  t-(m-i)  d/c  cos  <f>. 
When  the  optimum  spacing  is  used  in  calculating  t.,  we  denote 
the  t1  by  t..  1 


'Xt 

(5-58)  ^  -  ° 


Later  we  shall  see  that  these  relations  are  extremely  significant. 
From  the  general  expression  Eq.  (5-49),  we  also  obtain 


(5-59)  \  -  8„  =  ^  • 

w 


(T  *b<‘> 


vb  =  0 

which  are  independent  of  the  power  spectrum  of  the  interference. 
Notice  that  %  is  just  equal  to  the  first  term  of  Eq.  (5-49).  ^Since 
it  is  much  greater  than  the  second  term  of  the  same  equation,  «b  is 
approximately  equal  to  its  bound  . 

We  now  give  a  physical  interpretation  of  ab,  Bb  and  Yb*  The 
interpretation  will  also  apply  to  the  parameters  aa,  Ba  and  ya  for 
the  obvious  reason.  Under  hypothesis  H.  ,  the  output  of  the  sub¬ 
processor  Lbu(r)  may  be  rewritten  as 

’ 

Lbu(r)  =  j,  l)  [n1<V  +  wi(ti>]  •  ubi(ti>  •  dt 

-  JJ  T  ^.(t.)  +  w.(ti)j  •  ubj^zj)  *  *  dt  dz 

“b  Yb 

+  x  .  -£  -  y  ♦  -g- 

or  °r 

Thus  the  component  contributed  to  the  output  L.  (r)  by  the  desired 
signal  is  equal  to 

“b  Yb 

x  •  t  " y  •  t 

°r 

Similarly,  the  component  contributed  to  the  output  L.  (r)  by  the 
desired  signal  is  equal  to  D 


u  u 

y  *  t  ■ x  *  t 
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Apparently,  large  db  and  jlb  and  small  Tb  will  reduce  the  influence  of 
the  interference  in  making  decisions.  Referring  to  Eq.  (5-3),  the 
desired  signal  b(t)  may  be  decomposed  into  an  in-phase  component  and 
a  quadrature  component,  each  with  a  random  amplitude.  The  parameters 
dfc  and  jk  represent  the  bit  energies  contributed  to  Lbu(r)  and  Lbv(f) 
by  the  signals  Ub(t)  and  Vb(t),  respectively.  The  parameter  Tb 
represents  the  cross-energy  between  the  ub(t)  and  vb(t). 


Numerical  Results,  ASK  Signals 


In  this  section,  we  investigate  the  characteristics  of  the 
optimum  detector  for  detecting  the  on-off  keyed  binary  signals.  We 
first  investigate  how  the  system  responds  to  a  testing  signal  arriving 
from  an  arbitrary  angle  .  The  testing  signal  used  is  assumed  to 
have  the  same  waveform  as  the  desired  signal  b(t).  Again,  we  con¬ 
sider  the  case  when  the  interfering  signal  is  a  white  noise  process 
with  spectral  height  hn.  Later,  the  error  rate  performance  and 
methods  of  optimizing  this  performance  are  also  discussed. 


Letting  yk  ^  0  and  “b  %  &b  ™  E(l*  (5-34),  we  thus  obtain  the 
processor  for  detecting  the  on-off  keyed  signals  corrupted  by 
directional  white  noise: 


(5-60) 


Mr)  =  (<xfa  +  1)  Lbv(r)  +  U5  +  1)  Lbli(r) 


The  detection  system  makes  decisions  by  comparing  its  processing 
output  A(r)  to  a  given  threshold.  Whenever  the  processing  output 
is  less  than  the  system  signals  that  H  is  present;  otherwise, 
the  system  chooses  Hb. 

The  capability  of  a  detection  system  of  doing  its  job  is 
uniquely  decided  by  how  the  processing  output  under  one  hypothesis 
is  distinct  from  the  output  under  another,  and  the  value  chosen  for 
the  threshold.  A(r)  is  the  result  of  processing  the  interference- 
contaminated  inputs  in  the  manner  specified  by  the  right  side  of 
Eq.  (5-60).  Hence,  the  characteristics  of  this  detection  system 
are  determined  by  that  equation.  A(r)  is  randomly  distributed.  One 
important  characteristic  of  A(r)  is  its  mean,  which  we  now  examine. 

Since  the  te§ting  signal  b(t)  is  assumed  to  propagate  into  the 
array  from  angle  e,  the  input  to  the  ith  channel  of  the  processor 
is  given  by 
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(5-61)  =  n.(t.)  +  +  x  ubi(t.)  -  yvbi(t.) 

where 

(5-62)  ubl- (ti )  a  ub  Ct-(i-l)  ^  cos  e  -  (m-i)  ~  cos  <t>H 
and 

(5-63)  vbi (ti )  A  vbCt-(i-l )  ^  cos  0  -  (m-i)  ^  cos  <?] 


The  detection  system  is  optimized  for  the  desired  signal  arriving 
from  angle  e.  Hence,  the  sub-processors  Lbu(r)  and  l-bv(r)  will 
generate  the  outputs  as  follows: 


(5-64) 

a  ,  I  m  fT  . 

U.,(r)  =  JU  l 


-bu' 


Oi=S  J_T  ri(ti)ubi(ti)dt 


■5  f  fT  A 


’  uU  iJ-Tri(ti)uw(zJ)  ‘  f(t,z) 


and 


dt  dzl 


(5-65) 


(r)  = 


yi(t1)vb1<Vdt 


*'i(ti)vbj(zj)  *  ’ft*2)  ’  dt  dz 


Referring  to  Eq.  (5-61),  we  see  that  the  first  terms  of  Lbu(r)  and 
Lbv(r)  no  longer  represent  the^coherent  correlators,  because  the 
time-shiftings  in  ubi(t-j)  and  Vhi  (t-j )  are  not  equal  to  those  in 
ubi(ti)  and  Vbj(tj)»  respectively.  The  output  of  the  optimum 
processor  is  equal  to 


(5-66)  A(e,*;S)  =  (ab+l)tbzv(r)  +  (sb+l)  •  Cbu(r) 


Equation  (5-6 Q)  is  just  Jhe  special  case  of  A(e,<f>;e)  when  e=e.  The 
variances  of  Lbu(r)  and  Lbv(r)  may  be  determined  from  Eq.  (B-18) 
and  Eq.  (B-19),  respectively.  It  is  found  that 
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(5-67) 

var  Ubu(r)|Hb>  = 

1 

+ 

*2  ^ 
“b  + 

V 

and 

ar 

(5-68) 

var  <Lbv(r)|Hb}  - 

1 

7t 

(sb  + 

-2 

8b  + 

V 

where 

a2 

(5-69) 

V  V 

w 

2 

(5-70) 

A  A 

6b  —  IT"  <vb*  vb> 
w 

and 

2 

(5-71) 

A  .  0p  A 

Yb~  IT  <ub»  V 

Comparing  the  right  side  of  Eq.  (5-6)  with  that  of  Eq.  (4-2),  we  find 
that  they  are  the  same  if  $s=0.  We  also  find  that  v^(t)  is  equal 
{o  Eq.  (4-2)  if  ♦s  ■  -  tt/2.  Hence  the  closed  form  representation  for 
°*b  may  be  readily  obtained  from  Eq.  (4-11)  when  <t>  is  set  equal  to 
zero  and  s 

fT 

j  ea(t)  e(t)  .  dt 
is  changed  to 


As  a  result,  when  the  interference  is  a  white  noise  process  we  find 


(5-73)  3b  %  ab 

and 

(5-74)  Yb  %  0 

A 

Thus  when  the  interference  has  a  flat  spectrum,  the  mean  of  A(e,<j>;e) 
is  simply  given  by, 

(5-75)  E{A(e,(j>;0)  |Hb>  =  ^  (°5+1)  '  (°b  +  “b^ 

°r 

For  the  same  reason  as  given  in  Chapter  IV,  we  study  the  directional 
characteristics  of  the  system  by  examining  the  relative  angular 
response. 

E{A(e,<j>;0)  |H.  } 

(5-76)  F(e,4>;e)  =  20  •  log  — 

E{A(e,4>;6)  |Hb> 

The  signal  b(t)  is  a  random  process.  The  mean  of  its  energy  in 
the  observation  interval  [-T,TD  is  equal  to 


We  define  the  signal  energy  to  the  total  noise  power  density  ratio  as 
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is  the  signal  energy  to  the  internal  (white)  noise  power  density 
ratio.  We  proceed  to  investigate  the  directional  characteristics 
through  studying  the  angular  response  curve. 

Figure  5-3  shows  the  angular  response  of  the  detection  system 
(5°, 45°)  which  has  a  2-element  array  with  spacing  0.865a.  In  plotting 
this  curve,  we  assume  (S/N)  =  10,  h  |h  =  0.1  and 

W  il 

4  f  eb^^dt  = 

To  show  the  effect  of  varying  the  array  spacing,  the  response  of  the 
detector  (5°, 45°),  which  uses  the  same  number  of  elements  with  an 
increased^spacing  (1.73a),  is  given  in  fig.  5-4.  In  Fig.  5-4, 
F(5o,45°;0),.is  the  absolute  maximum  at  e=5°  and  <f>  becomes  a  relative 
minimum  at  0=45°.  Physically,  this  curve  implies  that  when  the 
signal  b(t)  propagates  into  the  array  from  the  angle  0,  the  average 
processor  output  is  the  highest.  On  the  other  hand  the  processor 
generates  almost  the  lowest  output  (in  average)  if  the  same  signal 
arrives  from  the  angle  <p.  As  a  matter  of  fact,  1.73a  is  the  optimum 
spacing  for  m=2,  0=5°  and  <t>=45°. 

In  Fig.  5-5  through  Fig.  5-8  we  show  the  angular  response  curves 
for  the  detectors  (55°, 45°)  which  are  designed  for  different  number 
of  antennas  and  different  spacings.  It  is  seen  in  these  curves  that 
the  response  at  the  angle  corresponding  to  the  incident  angle  for 
the  interference  is  always  lower  than  the  response  at  0=0.  This 
suggests  that  the  optimum  system  possesses  the  same  characteristics 
as  the  one  discussed  in  Chapter  IV.  In  fact,  F(0,<f>;0)  may  be 
written  simply  as 


vs.  0  given  m=2, 

Vhn  =  O'’- 


6  (  DEGREES  ) 


Fig.  5-7~F(55°,45°;e)  vs.  0  given  m=7,  d=0.537x, 
(S/N)=10  and  hw/hn  =  0.1. 


2 

(+  * 

— - ^ 

ab  +  ab 

The  parameter  ab  1S  often  much  greater  than  one,  hence  F(e,<f>;e) 
may  be  approximated  by 


(5-79)  F(0,4>;0)  $  40  •  log 


A  .  A  . 

Because  ab  and  mLa(0t4>;9)  have  the  same  directional  property,  it  is 
no  surprise  that  both  systems  have  the  similar  characteristics.  In 
plotting  Fig.  5-5  and  Fig.  5-7  an  optimum  spacing  is  used  by  each 
detector.  As  a  result,  the  average  processor  output  is  always  the 
highest  when  the  signal  b(t)  propagates  into  the  array  from  the 
desired  source. 


To  examine  the  behavior  of  F(0,<f>;0)  when  0=41,  we  show  four  curves 
in  Fig.  5-9  through  Fig.  5-12.  In  plotting  these  curves,  we  aqain  let 
(S/N)  =  10,  hjh_  =  0.1  and 


fT 


w' 


e?(t)dt 


1. 


)  (<J8) 


Fig.  5-12.  F(55°,5°;0)  vs.  6  given  m=4,  d=0.936x, 

(S/N)=10  and  hw/hn  *  0.1. 

A  A 

It  is  seen  that  F(e,<j>;0)  always  assumes  the  peak  value  at  0=e=<j>  for 
this  extreme  case.  When  m®#,  S/N=10  and 


the  parameter  ajj  is  equal  to  110  for  the  system  designed  for  the 
optimum  spacing.  For  the  case  shown  in  Fig.  5-9,  05  is  only  equal 
to  5.238.  The  detection  error  is  small  if  ab  is  large.  (The 
equation  for  detection  error  will  be  given  in  Eq.  (5-88).  Apparently, 
the  error  rate  performance  is  extremely  poor  when  e=<p. 

The  basic  properties  of  the  response  curve  discussed  in  this 
section  do  not  change  if  the  values  chosen  for  h#,  (S/N),  etc.  were 
different.  But  the  error  rate  performance  of  the  optimum  detector 
does  change  with  these  factors.  With  this  background,  we  next 
investigate  the  error  rate  performance. 

There  is  no  need  to  diagonalize  the  quadratic  equation  (5-34) 
when  the  interference  is  a  white  noise  process,  because  yb  £  0. 

But  nevertheless,  the  results  shown  in  Section  V.C  are  still 
applicable*.  Letting  yb  s  0  and  ab  *  Bb,  from  Eqs.  (5-38)  and  (5-39) 
we  obtain 


♦The  eigenvalues  are  Ci  ■  Btfiand  c2  =  “b  +  1  when  yb  -  0.  The 
new  and  the  primary  coordinates  are  related  by 


n 

f - 

O 

1 _ 

LbU<r> 

5bu. 

O 

_ 1 

LLbv<r>. 
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•■UV'V 


(5-80) 


2  2  “b  ,  .  1X2 

°bvb  0bub  ■  -?  (“b  +  '> 

o.. 


(5-81) 


a?  =  a?  =  Tr-(a,+l) 
bva  bua  T1  b  ' 

r 


Also,  as  a  result  of  Yb  £  0  and  %  Sb»  the  probabilities  of  making 
correct  and  false  decisions  on  whether  hypothesis  Hb  is  true  reduce 


(5-82)  Pd  $  e 


(5-83)  Pf  £  e 


In  practice,  the  probabilities  for  each  hypothesis  to  be  present  are 
equal  so  that  the  total  probability  of  making  wrong  decisions  on 

n  n 

— 7  — 7 

2af  2oh 

(5-84)  Pe  %  0.5  e  a  +  0.5  (  1  -  e  D  ) 

The  error  would  be  0.5  if  the  threshold  is  set  at  zero  or  at  infinity. 
Solving  3P~/3n  a  0  for  n,  we  find  that  the  lowest  total  detection 
error  by  adjusting  the  threshold  is  obtained  when 


(5-85) 


O,  aK 
o  a  b  i 

2  -7— ?  1 


V°a 


"(J) 


(<xb  +  1)‘ 


ln(ob+l)  a 
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J 
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At  n  =  we  denote  P<j  and  P^  by  Pd(r^)  and  Pf(no),  respectively. 

When  n  =  n  ,  the  two  probabilities  are  related  by 

O  \  Prf(n  ) 

(5-86)  Pf(V  -  4*  >><,< -  sV- 

a  b  b 

2 

where  <?a 

(5-87)  Pd(%)  -^jV<’b  -  (V')  “b 

The  corresponding  probability  of  making  wrong  decisions  on  H.  then  is 

__1_  D 

(5-88)  pe(n0)  =  0.5  -  0.5  •  ^rj  *  («b  +  D  “b 

b 

In  Eq.  (5-55),  the  parameter  a b  is  expressed  in  terms  of  x  and  v. 
The  delay  times  t  and  v  depend  on  the  incident  angles  e  and  <J>;  so 
does  the  optimum  spacing.  For  a  given  spacing  d  and  the  array  size 
m,  the  ratio  d/dop  varies  with  the  angles  e  and  4>.  Thus  the  parameter 
may  be  writtenKas  a  function  of  the  normalized  spacing  d/d. 


(5-89) 


n2 

_  %  r 

“b  "  V 


m  +  -pp- 

% 


s1nK"*j 


U  fT  ^ 


We  next  compare  the  detection  errors  Pe(np)  for  the  systems  use  the 
same  array  size  but  are  designed  to  be  optimum  for  the  different 
spacings.  We  shall  plot  the  error  pe(n«)  against  the  normalized 
spacing.  In  plotting  the  curve,  we  assign  values  to  (S/N)  and  hn/hw. 
By  doing  this,  the  signal  energy  to  the  internal  noise  power  density 
ratio  is  specified  alternatively  by 


(5-90) 


<(t) 


(S/N) 


■(■•a 


In  Fig.  5-13  we  show  Pe(no)  vs  d/dop,  where  0.2  <  d/dop  <_  2, 
for  m=2.  The  values  used  for  S/N  and  hn/hw  are  10  an?  31.61, 
respectively.  In  Fig.  5-13,  Pe(no)  becomes  minimum  when  the 
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Fig.  5-13 — Pe(n0)  vs.  d/dop  given  m»2,  hn/hw*31.61  and  (S/N)*10. 

normalized  spacing  Is  equal  to  1,  l.e.,  when  the  optimum  spacing  is 
used  by  the  detection  system.  The  same  type  of  curves  given  m=4  and 
8  are  shown  In  Fig.  5-14  and  Fig.  5-15,  respectively.  Again,  the 
detection  errors  P^r^)  assume  the  lowest  value  at  d  *  dQp. 

The  reason  for  the  detector  designed  for  an  optimum  spacing  to 
have  the  fewest  decision  mistakes  can  be  explained  mathematically 
by  the  behavior  of  Pe(no)*  In  Eq.  (5-88)  we  see  that  Pp(n0) 
decreases  as  ab  Increases.  Varying  the  spacing  d,  we  find  that 


d0p 


Fig.  5-14 — Pe(n0)  vs.  d/dop  given  m=4,  l^/hn-31.61  and  (S/N)=10. 

ab  reaches  the  maximum  at  dsd0p.  Hence  the  error  Pe(n  )  reaches  the 
minimum  at  this  particular  spacing.  0 

It  was  pointed  out  in  Section  V.D  that  ab  Is  proportional  to 
the  component  contributed  to  the  output  of  the  subprocessor  Lbu(r) 
by  the  desired  signal.  The  ability  for  the  Internal  noise  and 
external  Interference  to  affect  the  decisions  becomes  less  effective 
when  ab  gets  large.  Hence,  the  detection  error  is  minimized  when  an 
optimum  spacing  is  employed  to  design  the  detector.  For  the  rest  of 
this  section,  we  use  Pe(min)  to  denote  the  total  error  rate  for  the 
system  which  sets  the  threshold  at  r£  and  Is  designed  for  an  optimum 
spacing.  0 

To  emphasize,  Pe(mln)  Is  the  lowest  detection  error  we  can 
possibly  obtain  for  a  given  array  size  and  the  signal  energy  to  the 
internal  noise  power  density  ratio.  Moreover,  the  error  Pe(min)  is 
not  affected  by  the  incident  angles  e  and  4.  The  reason  is  that 
the  parameters  ab*  3b  and  yb  are  not  dependent  on  e  and  $  when  the 
system  is  designed  for  an  optimum  spacing. 
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Fig.  5-15— Pe(n0)  vs.  d/dQp  given  m=8,  hn/hw=31.61  and  (S/N)=10 


Up  to  this  point  In  this  section,  the  Interfering  signal  is 
assumed  to  have  a  flat  power  spectrum.  The  conditions  used  to 
get  Eq.  (5-56)  through  Eq.  (5-88)  are  Yb=0  and  “tr&b  which  are 
correct  because  the  Interference  Is  a  white  noise  process.  In 
Section  V.D  we  formally  show  that  regardless  of  the  power  spectrum 
these  conditions  are  automatically  satisfied  for  the  system  designed 
for  the  optimum  spacing  for  a  given  e,$  and  m.  Hence,  replacing 
“b  ^  (5-56)  through  Eq.  (5-88),  we  obtain  the  corresponding 

equations  valid  for  the  system  designed  for  the  optimum  spacing  and 
optimized  for  an  arbitrary  power  spectrum.  The  error  Pe(min) 
associated  with  this  detection  system  Is  equal  to 


(5-91) 


_  J_ 

(V1)  “b 


(5-92)  \  =  \  • 


Since  is  invariant  to  the  power  spectrum  of  the  interfering 
signal,  it  is  apparent  that  this  lowest  possible  error  Pe(min) 
does  not  change  from  one  power  spectrum  to  another.  Equation  (5-92) 
indicates  that  the  detection  error  Pe(min)  is  uniquely  determined  by 
the  array  size  m  and  the  signal  energy  to  the  internal  noise  power 
density  ratio. 

We  next  show  the  curves  of  Pe(min)  of  several  detectors.  In 
plotting  these  curves,  we  assume  that  hn  is  constant  so  that  the 
values  assigned  to  the  spectral  height  ratio  hn/hw  and  ratio  (S/N) 
specify  the  desired  signal  energy  to  the  internal  white  noise  power 
density  ratio.  For  a  given  hn/h*,  the  relation  between  (S/N)  and 
the  signal  energy  to  the  internal  noise  power  density  ratio  is 
given  by  Eq.  (5-90).  Since  we  shall  plot  Pe(min)  against  (S/N), 
the  relation  between  (S/N)  and  the  signal  energy  to  the  internal 
noise  power  density  ratio  is  shown  in  Fig.  5-16  for  convenience. 

The  error  Pe(min)  vs.  S/N  given  iip2  Is  shown  in  Fig.  5-17.  The 
values  chosen  for  the  spectral  height  ratio  hn/hw  are  3.161,  10, 

31.61  and  100.  For  a  given  S/N,  the  error  rate  Pe(min)  drops  as 
the  ratio  hn/hu  drops.  Physically,  this  means  that  the  lowest  possible 
detection  error  decreases  if  the  Internal  noise  level  decreases.  In 
Chapter  IV,  we  have  concluded  that  the  internal  noise  degrades  the 
performance  of  the  optimum  detector  L(R).  Here  we  find  that  the 
internal  noise  affects  the  detector  A(r)  In  the  same  manner.  For 
a  constant  internal  noise  level,  we  may  either  Increase  the  array 
size  or  raise  the  desired  signal  level  for  Improving  the  error  rate 
performance.  The  same  type  of  curves  given  m=4,  8  and  15  are  shown 
in  Fig.  5-18  through  Fig.  5-20.  It  Is  observed  that  P  (min)  drops 
steadily  as  the  array  size  increases  from  2  to  15.  e 

The  lowest  detection  error  available  to  the  detector  A(r)  is  far 
greater  than  the  error  provided  by  the  system  L(R)  for  detecting  the 
completely  known  binary  signals.  The  reason  for  having  such  a  large 
difference  is  obviously  attributed  to  the  random  phase  and  random 
amplitude  which  Inevitably  cause  many  decision  mistakes  for  the 
detect or  A(r). 


'V 
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Paulin)  =  0.5  -  0.5  •  — « 
e  &b  +  1 


where 


f  ek(t)*dt 
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CHAPTER  VI 

ON  SYSTEMS  DESIGNED  FOR  THE  OPTIMUM  SPACING 


A.  The  Interference-Free  Spacing 

When  the  desired  binary  signals  are  known  exactly,  we  have 
concluded  from  the  detection  error  given  in  Eq.  (3-32)  that  the 
number  of  wrong  decisions  is  minimized  when  the  detection  system 
is  designed  for  the  optimum  spacing.  The  same  conclusion  is  also 
obtained  for  the  random  phase  and  amplitude  ASK  binary  detector 
discussed  in  the  last  chapter.  In  this  section,  we  examine  L(R), 
Lbu(f)  and  Lbv(0  for  a  unified  explanation  for  this  phenomenon. 
It  will  become  clear  later,  the  optimum  spacing  may  be  called  the 
interference-free  spacing. 


The  definitions  for  the  random  variables  L(R),  Lbu(r)  and 
Lbv(r)  are  very  similar.  They  may  be  represented  by  the  general 
expression  L  given  below 


(6-1) 


m 

l  r1(t1)g1(t1)  •  dt 
i=l  1111 


-  I 

t- 1 


1  1 


m  + 
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-T 


♦,<*> 


A  _ — h 
A  W 


m 

l  r . (t. )dt  . 

i=l  1  1 


rT  m 

I  94(2,-)  *  dz 

-T  1  j=l  3  3 


because  L(R)  =  Lg  if  g(t)  =  s(t),  Lbu(r)  =  Lg  if  g(t)  =  Ub(t)  and 
Lbv(r)  '  Lg  when  g(t)  =  -Vb(t).  When  optimum-spaced  arrays  are  used 
in  both  detection  systems,  we  have  shown  in  Section  IV. F  and  Section 


V.D  that 

m 

g1(%i)  -  0 

(6-2) 

A 

and 

(6-3) 

m 

,1, 

ubi(ti>  vbl(ti>  =  0 

1 
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In  Eq.  (6-1),  as  a  result  of  using  an  optimum-spaced  array,  the 
second  term,  which  reflects  the  characteristics  of  the  interference, 
vanishes.  Hence, 


m  rT  m 

(6-4)  L(R)  =  l  r(t-)  s-(t-)  •  dt 
J_T  i=l  1  1  11 


(6-5)  £^(r)  •  h^1  f  r,(t,)  u^)  •  dt 

and 

(6-6)  tbv(D  -  -h;1  [T  .1  r,^)  vbim  .  dt 

Equation  (6-4)  through  Eq.  (6-6)  imply  that  the  sub-system 
which  varies  with  the  power  spectrum  of  the  interfering  signal  is 
no  longer  needed  when  the  detection  system  is  designed  for  the 
optimum  spacing.  The  input  to  the  ith  channel  under  hypothesis  Hb 
is 


(6-7)  rftj  =  n  Ct-(m-l)v]  +  w^)  +  b^fc.) 

where  the  signal  b(t)  is  assumed  to  be  deterministic  if  rf(fcj)  is 
substituted  into  Eq.  (6-4),  and  b(t)  is  assumed  to  have  a  random 
phase  and  amplitude  when  rj(tf)  is  substituted  into  Eq.  (6-5)  and 
Eq.  (6-6). 

(6-4)  and  simplifying  the  equation 


fT  m  i\,  % 

.  dt  +  l  b-(t.)  w^t.)  •  dt 

J-T  1=1  1  1  11 


Similarly, 


Substituting  r-f (%-j )  into  Eq. 
by  using  Eq.  (6-2),  we  obtain 

(6-8) 

i\j  fT  m  ^ 

*-(R)  *  I_t  -li  bi^i^  si ^i 


and 


(6-10) 

4,„<r>  -  y 


Vbi(fyWi(fy  ’  dt 
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The  simplified  versions  of  L(R),  Lbu(r)  and  Lbv( r)  do  not 
contain  the  external  interfering  signal  n(t).  This  implies  that  the 
influence  from  the  interference  in  the  decision  process  is  eliminated 
completely  by  the  processor  when  the  detecto^  is  designed  for  an 
optimum  spacing.  When  the  processor  output  L(R)  is  compared  to  the 
threshold  to  make  decisions,  only  the  internal  white  noise  w-j(t)  can 
cause  decision  mistakes  for  the  case  when  the  desired  signals  are 
completely  known.  For  signals  with  random  amplitude  and  phase,  the 
uncertainty  in  the  values  of  x  and  y  also  cause  decision  errors.  This 
explains  why  the  detection  error  for  the  optimum  system  a (r)  is  far 
greater  than  the  system  L(R).  But,  nevertheless,  all  of  l(R),  ^(r) 
and  Lby(r)  are  interference  free.  Physically,  we  may  say  that  the 
optimum  spacing  is  an  interference-free  spacing.  For  a  given  m,  hw, 
etc.,  it  is  this  very  property  which  enables  the  detector  designed 
for  an  optimum  spacing  to  have  the  lowest  detection  ej^ror.  Because 
the  eigenfunctions  do  not  appear  in  L(R),  Lbu(r)  and  Lbv(r)  any  more, 
the  structures  for  detectors  L(R)  and  fr(r)  do  not  vary  from  one  power 
spectrum  to  another. 

B.  FSK  Binary  Detector 


Equation  (5-2)  represents  FSK  binary  signals  when  u^b*  In 
practice,  the  carrier  frequency  under  one  hypothesis  differs  from 
that  under  the  other  only  slightly.  Also,  it  is  a  general  practice 
to  use  identical  modulating  signals  for  both  hypotheses.  Let 


(6-11)  ek(t)  *  e(t),  k  =  a  and  b, 
and 

(6-12)  <*>k  -  «0  +  Vd  . 

where  a>d  «  up,  ej<  =  -1  under  hypothesis  H  and  ek  *  1  under  hypothesis 
Hb.  For  simplicity,  we  define  a  K 

(6-13)  u0(t)  A.  e(t)  •  cos  u)Qt 


T 


(6-14)  v0(t)  A  e(t)  •  sin  u>Qt 

Then,  from  Eq.  (5-6),  we  have 


(6-15)  uk(t)  *  uQ(t)  *  cos  “d*  "  e|(  *  V0U)  *  sin  “dt 
and 


(6-16)  vk(t)  *  ek  *  uQ(t)  *  s1n  “d*  +  v0^  *  cos  “d1* 


In  this  section,  we  Investigate  the  form  of  the  optimum  system  for 
detecting  the  FSK  binary  signals  in  Eq.  (5-2)  when  the  spacing  chosen 
for  the  array  Is  equal  to  the  optimum  spacing  defined  In  terms  of  the 
wavelength  A  associated  with  the  angular  frequency  x.  More 
precisely,  0  0 


(6-17)  d  A - ^ - 

p  m  |cos  e  -  cos  <fr| 


At  this  particular  spacing,  if  jcos  e  -  cos  <|>|  is  not  too  small, 
negligible  error  will  be  Introduced  If  we  approximate 


(6-18)  u^Ct  -  (m-l)v  -  (1-1)(t-v)D  %  u^t 

A/  <\f 

Consequently,  uk^(tj)  and  vk^(t,j)  may  be  approximated  by 

(6-19)  Uk1  (t-j )  ^  uoi(*i)  *  cos  “d*  "  ek  *  voi^1^  *  sin  “d* 
and 


(6-20)  v^ftj)  £  ek  •  uQi(ti)  •  sin  «dt  +  vQi 

m  <v  m  <v, 

respectively.  Since  l  s  l  voi^i^  * 

from  Eqs.  (6-19)  and  (6-20)  that 


A 

(t^)  •  COS  Wjt 

0,  we  conclude 


(6-21) 
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uki*V 

*  0 

m 

(6-22) 

vki(ti> 

=  0 

and 

(6-23) 
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A# 

*v 

A 

uki*V 

vki<t1>  =  0 

The  eigenfunction-dependent  terms  once  again  drop  out  of  the 
equation  of  the  test  variables  Lku(r)  and  Lkv(r);  so  we  find 

*\i  *i  fT  m  a 

(6-24)  Lku(r)  =  h;1  ^  r^t,)  u^tt,)  •  dt 

and 


(6-25)  Lkv(r)  =  -h;1  ji  rl(t,)  v^t,)  •  dt 

where  u|<i(tf)  apd  vkj(t^)  £re  given  by  Eq.  (6-19)  and  Eq.  (6-20), 
respectively.  Lku(r)  and  Lkv(r)  can  be  shown  free  from  interfering 
signal.  For  this  reason,  it  is  believed  that  the  FSK  binary  detector 
designed  for  an  optimum  spacing  also  provides  the  lowest  detection 
error.  Because  an  identical  modulating  signal  is  employed  under 
each  hypothesis,  we  know  that 


(6-26) 


-v-  “v  r 

“k  =  *k  s  TT 

Wf 


k  *  a  and  b. 


Hence, 


(6-27) 


Ua(a.0»Y)  *  Mb(a»e»Y) 


Resetting  the  threshold,  from  Eq.  (5-24),  we  obtain  a  reduced 
likelihood  ratio  test  for  FSK  binary  signals  received  by  an 
optimum-spacing  array. 


□ 
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(6-28)  iLk(r)  A  Lku(r)  *  tL(r)  -  £,(r>  -  L *,<r) 


<  threshold 

H. 


v  %  'V 

Writing  Lku(r)  and  L|<v{r)  directly  in  terms  of  u0ii[ti)  and  v0i(ti), 
and  substituting  the  new  equations  for  tku(r)  and  L^yCr)  into 
Eq.  (6-28),  we  obtain  a  simpler  optimum  detector: 


(6-29) 


fT  fT  r\j 

Afsk^r)  =  VV  *  sin  “d*  *  dt  *  W  *  cos  U)dt  *  dt 

■'-T  ;-T 

T  %  tT  ^ 

-  j  Au(t.)  *  cos  •  dt  •  j  \(t^)  •  sin  u)dt  •  dt 


where 


<  threshold 
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f\j  111  r\, 

W  ’  ,1,  ri(V  uoi(ti> 


%  m  % 


(6-30)  Ay(ti)  =  J  r^t^  vQi(t1) 


B.  On  General  Random  Phase  and 


In  this  section,  we  shall  briefly  discuss  the  optimum  detector 
for  detecting  the  binary  signals  whose  amplitudes  and  phases  are 
distributed  differently  from  what  were  assumed  In  the  last  chapter. 

In  theory,  the  joint  probability  density  function  of  the  random 
variables  x  and  y  defined  In  Eq.  (5-4)  and  Eq.  (5-5)  may  be  de¬ 
termined  once  the  statistics  for  the  amplitude  and  the  phase  are 
specified.  Let  the  joint  density  be  p(x,y).  The  density  of  R  is 
obtained  by  averaging  the  conditional  density  Eq.  (5-12)  over  x 
and  y.  For  simplicity,  we  define 


106 


-x2aj-y26£+2xyy£ 

jjp(x.y)  •  e  +  x  L ku(r)  +  y  Lky(r)  •  dx  dy 

A  pCLku(r).  Lky(r),  “t»  8k*  Yk*»  cXy^ 

where  aj^  3k  and  yk  are  equal  to  ak*  ek  and  respectively  If 
a2  =  1,  ana  CXy  represents  all  parameters  which  characterize  the 
random  variables  x  and  y.  Then  the  density  of  R  under  hypothesis 
Ht,  k=a  and  b,  is  equal  to 

Lr(r) 

p(R|Hk)  *  lim  —===■  -paku(r),  Lkv(r),  sj,  v{,  C  ] 

N~  J(2.)mN  |K| 


and  the  likelihood  ratio  test  takes  the  following  general  expression 


(6-31) 


LRT  = 


_  Pabu(r)*Lbv(r)»  °£»  •••«  CXy^ 


Paau<r)’Lav(r>’  a'a  *’  V 


-fabu^r^Lbv^r)**“*Lau(r^Lav(r^***ab»aa»  ***'  V 
<  threshold 

If  we  implement  the  function  fCO,  we  obtain  the  optimum  detection 
system.  Since  fCO  is  expressed  as  a  function  of  a£,  Lhu(r),  Lav(r), 
etc.,  evidently,  the  sub-processors  which  compute  Uu(r)  and  Lkvlr), 
k=a  and  b,  for  the  system  Ag (r)  are  again  required  by  this  general 
optimum  detection  system  Eq.  (6-31).  In  other  words,  the  basic 
building  blocks  for  the  random  phase  and  random  amplitude  binary 
detector  do  not  change  with  the  assumptions  on  the  amplitude  and 
phase.  But  how  the  outputs  of  these  sub-processors  should  be 
processed  is  determined  by  the  function  p(x,y)  because  fC*3  does 
vary  with  p(x,y).  We  have  demonstrated  that  the  variables  Lku(r) 
and  Lkv(r)  are  free  from  the  interference  if  the  detector  A(r)  is 
designed  for  an  optimum  spacing.  Even  though  the  formula  for 
the  detection  error  of  the  system  fCO  Is  not  known,  there  is  a 
reason  to  believe  that  for  a  given  array  size  m,  spectral  height 
h#,  etc.  the  lowest  detection  error  will  be  provided  by  the  system 
designed  for  an  optimum  spacing. 


107 


If  we  set  <t>s=0  in  Eq.  (5-2),  then  the  binary  signals  would 
contain  one  random  parameter  which  is  the  amplitude.  From  Eqs. 
(5-4)  and  (5-5)  we  find  that  y=0  when  *  *0.  Hence  the  conditional 
density  for  R  under  Hk  becomes 


(6-32)  p(  R|x,Hk) 


Lr(r) 


2  ' 

x  d 


lim 

N-~° 


+  xLku(r) 


J  (2tt  ) 


mN 
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Accordingly,  the  likelihood  ratio  test  associated  with  the  random 
amplitude  binary  signals  is  equal  to 


2  . 


X  a 


(6-33>  fp(x)  .  e'  *  X  Lb“<r)  •  * 


LRT  = 


2, 


-X  a 


|p(x) 


+  2  Lau(r)  •  dx 


Hence,  the  sub-processors  which  compute  Lau(r)  and  Lbu(r)  in  the 
system  Ag(r)  are  the  building  blocks  for  the  random-amplitude 
binary  detector. 


For  most  cases  of  interest,  the  analytic  solutions  for  the 
likelihood  ratio  test  are  not  obtainable.  As  a  result,  various 
types  of  approximations  may  be  needed. 
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CHAPTER  VII 
CONCLUSIONS 


The  goal  of  this  research  was  to  find  the  optimum  detection 
system  for  detecting  corrupted  binary  signals  received  by  a  linear 
antenna  array  in  an  effort  to  minimize  the  detection  error.  The 
binary  signals  may  either  be  completely  known  or  have  a  random 
amplitude,  which  is  Rayleigh  distributed,  and  a  random  phase, 
uniformly  distributed  between  -it  and  u.  The  approach  used  to 
determine  the  optimum  system  was  to  perform  the  likelihood  ratio 
test  based  on  the  statistical  detection  theory. 

The  results  presented  here  are  obtained  from  analytic  solutions. 
It  has  been  possible  to  obtain  these  solutions  mainly  because  we  have 
derived  an  analytic  expression  for  the  inverse  of  the  covariance 
matrix  K.  Hence  this  research  removes  a  mathematical  difficulty 
encountered  frequently  in  analyzing  problems  related  to  optimum 
space-time  signal  processing. 

In  the  completely  known  signal  case,  the  desired  signals  con¬ 
sidered  in  Chapter  III  are  the  most  general  ones  in  the  sense  that 
the  waveforms  may  have  arbitrary  expressions.  Hence,  the  equations 
formulated  for  implementing  the  optimum  processor  and  for  de¬ 
termining  the  system  performance  are  valid  for  the  ASK,  FSK  and 
PSK  binary  signals. 

For  the  case  when  the  amplitude  and  phase  of  the  desired  signal 
are  random,  the  equation  needed  for  Implementing  the  optimum 
processor  was  also  derived.  Since  the  phase  was  assumed  uniformly 
distributed  between  -it  and  tt,  the  result  can  not  be  applied  to  the 
PSK  binary  signals.  In  general,  we  shall  have  to  use  a  numerical 
method  to  determine  the  error  rate  performance  of  the  system  when 
the  signal  has  a  random  phase  and  amplitude.  But  when  the  binary 
signal  is  an  on-off  keyed  type,  a  simple  analytic  formula  for  the 
detection  error  Is  available. 

The  optimum  processors  for  both  cases  are  basically  cor¬ 
relation  receivers.  In  general,  the  structures  for  both  types 
of  processors  are  dependent  on  the  power  spectrum  of  the  Inter¬ 
fering  signal. 

The  directional  characteristics  of  both  detection  systems  have 
been  investigated  through  examining  how  the  systems  respond  to  a 
testing  signal  arriving  from  an  arbitrary  angle.  We  have  found 
that  both  processors  respond  unfavorably  to  the  testing  signal  in 
the  sense  that  the  processor  outputs  are  suppressed  If  the  signal 
propagates  Into  the  antenna  array  from  the  angle  4>.  When  the 
testing  signal  arrives  from  the  desired  angle  e,  the  processor 
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outputs  are  high.  The  ability  for  each  processor  to  favor  the 
signal  from  the  angle  d  and  to  suppress  the  signal  from  the  angle 
depends  on  the  angular  separation  between  e  and  <t>.  It  has  been 
shown  that  this  ability  can  be  Improved  by  either  increasing  the 
array  size  m  or  by  adjusting  the  spacing  d  if  the  angles  e  and  Vs 
are  fixed. 

The  detection  error  is  determined  by  the  array  size,  the 
internal  noise  level,  the  desired  signal  power  level  and  the 
spacing  d  for  which  the  detection  system  is  designed  to  be 
optimum.  Vie  have  concluded  that  the  detection  error  is  minimized 
If  the  array  elements  are  spaced  at  an  optimum  spacing,  which  is 
solely  determined  by  0,  $  and  m.  The  reason  for  this  is  that 
the  processor  outputs  are  interference-free  when  this  particular 
spacing  is  chosen.  Another  advantage  for  using  an  optimum-spaced 
array  is  that  the  power-spectrum-dependent  processing  units  can 
be  removed  from  the  processor  without  affecting  the  error  rate 
performance.  Hence,  for  a  given  array  size,  signal  power  level, 
and  the  spectral  height  h#,  the  optimum  detection  system  is  also 
Invariant  to  the  power  spectrum  of  the  interfering  signal. 

In  a  realistic  situation,  there  are  always  at  least  some 
parameters  characterizing  a  communication  signal  that  are  not 
known.  In  Chapter  III  we  have  idealized  the  situation  by 
assuming  the  binary  signals  to  be  completely  known.  The  detector 
thus  obtained  represents  an  ideal  case.  The  error  rate  associated 
with  this  case  will  be  lower  than  that  for  signals  with  random 
parameters.  Furthermore,  this  ideal  detector  is  optimum  in  the 
sense  that  for  a  given  condition  (e.g.,  array  size,  spacing  d, 
etc.)  its  detection  error  represents  a  lower  bound  for  the 
performance  that  might  be  obtained  with  other  types  of  array 
processing  (such  as  with  adaptive  array  techniques  C24-26]). 

That  is,  the  error  rate  given  by  Eq.  (4-30)  represents  a  lower 
bound  for  the  detection  error  for  all  conceivable  binary  detectors. 

When  the  phase  and  amplitude  are  distributed  otherwise,  we  have 
formally  shown  that  the  new  likelihood  ratio  Is  expressible  in  terms 
of  random  variables  L|<u(r)  and  Lkv(r)  and  the  parameters  a^,  Sk  and 
Yk»  For  this  reason,  discussions  given  in  the  latter  portion  of  this 
research  are  believed  quite  general. 
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APPENDIX  A 

A  PROOF  OF  EQ.  (2-36) 


In  this  appendix,  a  simple  proof  of  Eq.  (2-36)  is  given.  K*  and 
Q*  are  the  covariance  matrix  and  its  inverse,  respectively. 


( A— 1 ) 


(A-2) 


K*1 

K*2 

•"  K*m 

K*  * 

K|! 

K22 

"•  “Si 

• 

• 

• 

• 

• 

• 

• 

• 

• 

JBi 

•Sz 

...  K* 
mm. 

Ifl 

<Vf 

Qlm 

«!t 

Ofz 

—  Q2m 

Q*  * 

• 

• 

• 

• 

• 

• 

• 

• 

• 

<Si 

to 

•  •  •  0* 
inm 

The  general  expression  for  the  submatrix  G*.  of  the  product  of  K* 
and  Q*  Is  equal  to  1J 

(A'3)  k*p  qpj-  1  4  '•  j  s  ■ 

Me  need  to  prove  that  G*.,  is  the  null  matrix  if  i f*= j  and  G*.  is  the 
identity  matrix.  It  has J been  pointed  out  that  both  and  QT.  are 
diagonal,  i.e. ,  j  j 


The  product  of  K?.  and  Q?.  is  also  diagonal.  In  terms  of  the  elements 
of  K*.  and  Q?.,  Jg?i..,  Jthe  element  of  GT.  at  u  position  is  equal 

£q  ijxx  ij 

m 

(A"6)  g1jns  kf(i-l)N+£X(p-l)N+)i]  qC(p-l )N+£X(j-l  )N+£l 

From  Eqs.  (2-28)  and  (2-33),  we  have 

(A-7)  kf(i-l)N+OC(p-l)N+il]  =  XS.  '  (1_5ip)hw 

Substituting  Eqs.  (A-7)  and  (2-36)  into  Eq.  (A-6),  yields 

m  /  .  /  5 . . 

(A-8)  ^,(X*  -(,-4ip)hw)  ‘(tJ  ' 

For  i^j,  we  have 

g,  .  -("'-zMW,)2  MW 

+  ^£“hw^X£+^m_Z^X£“hW^  q 

hw^r^-Dhw3 

Hence,  Gt^,  i^j,  is  equal  to  the  null  matrix.  If  i=j,  we  have 

.  -<m-’)(xrhw)2  _  »tC  xt*(m-2)(»rhw)] 

i  0  0  _ 

hwCmA 4- (m- 1 ) hw3  hwQnA  z- (m- 1 ) hw] 

=  1  . 

Hence,  G*.  equals  the  identity  matrix.  We  complete  the  proof  that 
Eq.  (2-3o)  is  correct. 


_ \ 

hwDnX1-(m-l)hwV 
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APPENDIX  B 

THE  STATISTICS  OF  L(R),  Lku(r)  and  l_kv(r) 


In  this  appendix,  we  derive  the  variances  and  formulate  some 
useful  relations  concerning  L(R),  Lku (r)  and  L|<v(r)  defined  in 
Chapter  III  and  Chapter  V.  The  expressions  for  them  are  very  similar. 
In  order  to  avoid  repetitions,  we  define  two  new  functions  and  Lg: 

(B-l ) 


(B-2) 

Yi(t)  = 


f  n^ (t)  +  w.(t)  +e  a^(t),  under  Ha 
n^(t'  +  w..(t)  +e  b^(t),  under  Hb 


The  signals  a(t)  and  b(t)  may  contain  random  parameter  r.  Because 
the  equations  to  be  derived  will  not  be  directly  expressed  in 
terms  of  a(t)  and  b(t),  we  place  e,  e=0  or  1,  in  Eq.  (B-2) 
to  see  the  effects  when  a(t)  or  b(t)  equals  zero. 

For  a  given  r,  the  means  of  and  under  Hb  are 
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“  I  I  b.(t.)  •  J  J  I*  v(t,z)  •  dt  dz 

JLt  hM 

The  product  of  two  conditional  means  is  equal  to 
(B-4) 


E{Lf|r,Hb)  •  E(L  |r,Hb) 


m  rrT 


-  c 


.y,  JJ_T  W  •  W  •  W 


dt  dz 


m 


.  -i 


bjL(x£)bi(ti)-f£(x1)gj(zj)*'i'(t,z).dx  dt 


dz 


m 


-«*  is, 


rT 


biL(xx,)bi (ti)*gJL(xSL)fj(zj)*li'(t»z)'dx  dt  dz 


m 


+  c  .  mi 


_T  WMV  •  V9-,’  •  9j(zj>  • 


H»(x,y)  •  4>(t,z)  •  dx  dy  dt  dz 


The  product  of  Lf  and  L  is  given  by 


(B-5) 


LfLg 


m  rr  I 

=  ll  Y^tJ  Yj(Zj)  •  f,(t.)  g.(z.)  •  dtdz 

i,j=l  JJ-T  1  1  3  3  1  1  3  3 

m  rrrT 

-  m  Y£(xJ,)Yi(ti)‘fj,(x  )g.(z.)*y(t,z)-dx.  dt  dz 

i,£,  n=VJJ-T  *  *  1  1  JO  X.  j  j 


m  rrrT 

-  m 

£,i,j=l  J 


;; 


T^(x£)Yi(ti)*fj(zj)g,(x£)*'1'(t»z)-dx-dt  dz 


m 


•  v(x,y)  •  dx  dy  dt  dz 

For  a  given  r,  the  mean  of  yAtA)ya(xa)  under  hypothesis  H.  is  given 
by  i  i  £  £  b 

(B— 6)  YA(x£)|r,  Hb} 

=  R(t-x)  +  hw  .  «l4  •  fi(t-x)  +  e2  .  b1-{ti)  b£(x£) 

Hence,  the  conditional  mean  of  L^L  under  hypothesis  Hb  is  given  by 
(B— 7 ) 

E{Lf  Lg|r  ,Hb} 

m  rrl 

*  |J.T  R(t-2)  •  9j(2j>  •  dt  dz 

'  "  A,  £  R(t'*)-CV\>V2j)+Vxii)V24)]  • 


^(t.z)  •  dx  dt  dz 


(B-7)  +  af  HI  R(t-x)  •  fAz.)  g  (y  ) 

(Cont)  j,n=l  ""-T  3  3  n  n 


y(t,z)  •  y(x,y)  •  dt  dz  •  dx  dy 
m  t  T 


m  I 

+  "w  ,1,  J.T  W  91(tt)  •  dt 

m  r  tT 

-  h  11  Cf - (t . )g -(z .)  +  g.(t.)f ,(z.)]»’i'(t,z).dt  dz 

w  i,j=l  1  1  J  J  i  t  J  J 

m  rrrT 

+  n,h  *  ll  (y  ).y(t,z).y(t,y)*dt  dz  dy 

w  j,n=V"-T  3  3  n  r> 

9  m  f  fT 

+  ‘  1  Jl  JJ-T  b<(t1)b3(ZJ)  '  f1  (t1 }  Sd(2J>  ’  dt  dZ 

•  y(t,z)  •  dx  dt  dz 


Slj=l  JJJJ-T^1^0^  *  fn(yn)9j(Zj}  * 


y(x,y)  •  y(t,z)  .  dx  dy  dt  dz 


The  expression  given  above  may  be  simplified  by  applying  the 
following  relations: 


(B-9)  R(t-z)  •  l_  (it-hw)  •♦eCt)*t<z) 

The  second  term  of  Eq.  (B-7)  Is  equal  to 
(B-10) 

m 

-2m  •  H  R(t-x)f  (xf)g.(z.)  .  'f(t.z)  •  dx  dt  dz 
j,*=l  *  1  3  J 


m  r  rT  oo 

-  -2m  II  I 

i,j=1  "-T  1  1  J  0  1=] 


VhW 

- f c - dz 


m  + 


w 


x£  “  hw 


The  third  term  of  Eq.  (B-8)  equals 
(B-ll ) 


l!  jf  g^)^)  *  [- 

,j=l  JJ-T  1  1  J  J  *=1 


m 


f.j 


x  -  h 

h  "  dz 


The  fifth  term  of  Eq.  (B-8)  equals 
m  rrT 


.  m  rrl 

(B— 12)  -2  h^lj  ^  Jj  ^  f^tj)  9j(Zj)  •  <Kt,z)  •  dt  dz 


The  sixth  term  of  Eq.  (B-8)  is  equal  to 
(B-13)  m 


vv  ffT  “  ♦  .( t)4>  (z) 

• h- &,  JJ.T  WW-  K  \2 dt  d2 

l  v^V 


Substituting  Eqs.  (B-10),  (B-ll),  (B-12),  and  (B-13)  into  Eq.  (B-7) 
and  simplifying  the  resultant  expression,  yields 


(B-14) 


E(LfL  |r.  Hb) 


V  j,  fT  W  91(t1)dt  '  ,g.,  |fT  fi(t1)9j(2j>  • 

k 


•  Y(t »z )  •  dt  dz 


,  J,  |fT  •  W  W  • dt  dz 

HI.  }|f  ",(»,)•>,(*!)  •Cf1(xt)gj(zJ)+gi(xt)fj(2j)] 


2 

-e  • 


•  v(t,z)  .  dx  dt  dz 


+  e  M  ,,  JJJJ  T  W  MX£>  •  W  9j(zj)  * 


•  ^(x,y)  •  ^(t.z)  .  dx  dy  dt  dz 

The  sum  of  the  last  three  terms  In  Eq.  (B-14)  are  equal  to  Eq.  (B-4). 
Case  1 

We  next  determine  the  variance  of  L(R).  If  and  Lg  are  equal  to 
L(R)  if  f(t)  *  g(t)  »  s(t).  Since  the  signal  b(t)  does  not  contain 
random  parameters,  the  conditional  mean  is  just  the  mean.  Hence  the 
variance  of  L(R)  under  Hb  is 

(B-15) 

var{L(R)|Hb}=  E  {L2(R) |Hb)  -  E2{L(R) (Hb> 

"  {Eq.  (B-14)}f(t)=g(t)=s(t)  -  {Eq.  (B-4) >f(t)*g(t)*s(t) 
m  »T  n  m  r  r T 

"  VJ,  |.Tsf(ti>dt  -  ||.T  dz  • 

Ik 
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The  variance  of  L(R)  under  Ha  may  be  derived  In  the  same  manner. 

We  have  found  that  the  variance  under  H.  has  the  same  expression. 

O 

Case  2 

When  b(t)  *  x  ub(t)  -  y  vK(t),  where  x,  y,  uk(t)  and  vk(t)  are 
defined  In  Chapter  V,  then  D  do 


(B-16) 


E{bi(ti)  Vzj»  a  ^Cubi^i)  W  +  vbi(V  vbo(2j)I] 


E{b.(t.)}  =  0 
From  Eq.  (B-3),  we  obtain 


E(Lf lHb>  ■  E{E{Lf|r,  HbH  *  0 


Similarly, 


HL  |Hb)  =  0 


For  the  sake  of  simplicity,  we  define  that 
m  fT  m  ffT 


in  ft  m  rt  I 

».f>  -  J  T  ui(ti)fi(ti)dt  -  ll  ^  jj  ^  ui(ti)fj(z:j)4'(t,z).dt  dz 


in  1 1  m  , 

<v,g>  a  l  vi(ti)g.(ti)dt-  JT  v.(t.)g.(z.)y(t,z)-dt  dz, 
1=1  J-T  1  1  1  1  175*1  JJ. j  1  1  J  J 


etc.  From  Eq.  (B-14)  and  Eq.  (B-16),  we  obtain 
(B— 17) 

E<Lf  LglV  *  E{E{Lf  Lglr»  Hb» 

2  2 

3  \  •  <f»9>  +e  or  <ub*f>  •  <ub,g> 

2  2 

+  e  o„  <vK,f>  •  <vk,g> 


i 


E{Lf*Lg|Ha}  is  readily  obtained  when  the  subscript  b  in  Eq.  (B-17)  is 
changeu  to  a.  The  parameters  ak,  3k  and  iy,  defined  in  Chapter  V  may 
be  written  as 


“k  =  IT  <uk*  uk>  *  k  =  a  and  b 

nw 


8k =  i r  <vk»  V 

w 

o2  o2 

l  ’  ^  <uk>  V  *  <vk*  V 


1 

1 

*1 


1 


3 


We  next  formulate  the  variances.  Lf  and  Lg  are  equal  to  Lku(y)» 
k  =  a  and  b,  if  f(t)  *  g(t)  ■  h“‘  uk(t).  Hence, 

(B— 18) 

var{Lku(Y)|Hk}  =  E{L2u(Y)|Hk)  -  E2{Lku(Y)  |Hk> 


<uk'V 


+ 


e2°?*hw2<uk'uk>2  +c2°?h;2<uk»vk>2 


.  2,  2 
+  e  (ak 


Lf (y)  and  L  (y)  are  equal  to  L.  (y)  if  f(t)  =  g(t)  =  -h'V (t). 
Hence  g  KV  w 


(B-19)  var{Lkv(Y) |Hk>  =  ^f{8k  +  e^(3k  +y  k)> 

To  examine  if  Lku(Y)  and  *-kv?Y)  are  uncorrelated,  we  compute 
(B— 20) 

E{Lku(r)Lkv(Y)lHk>  -  E{Lku(Y)|Hk}  E{Lkv(Y) |HkJ 


s  ’  T  {Yk  +  ^kK  +  6k)} 
ar 

Hence,  in  general,  L)cu (y)  and  L^y)  under  H|<  are  correlated.  If  an 
optimum-spaced  array  is  used  in  the  detection  system,  we  have  shown 


'•L. 
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in  Chapter  V  that  sJ).  When  the  interference  is  a  white  noise 
process,  we  know  that  y^  0.  From  Eq.  (B-20),  we  conclude  that 
L|<u(y)  and  Lj^y)  are  uncorrelated  for  those  two  special  cases. 
Other  relations  can  also  be  obtained  from  the  general  formula 
Eq.  (B-17).  For  example, 

(B-21 ) 

E{Lau(Y>Lbu<Y>lHb} 


-1 


e2o? 


=  hw  <ub,ua>  +  -J-  <u5.ub><ub,ua>  + 


2  2 

e  ar 


^b’V'W 


w 


a’V 


+  e  a. 


<ua,ub>  +  e-y, 


b<ua’V 


Hence,  Lau(y)  and  Lby(y)  are  also  correlated.  In  Section  V.B,  we 
need  to  know  the  variances  of  Lau(y)  and  Lav(y)  when  e=  0.  Letting 
c=0  in  Eq.  (B-18)  and  Eq.  (B-19J,  we  obtain  the  variances: 

(B— 22)  var{Lau(y)|Hk}  =  ^ 


(B— 23)  var{Lav(y)|Hk)  * 

°r 

We  can  show  from  Eq.  (B-17)  that  L|cu(y)  and  Lqy(y)>  where  k,  p 
and  b,  are  uncorrelated  if  an  optimum-spaced  array  is  used  in  the 
detection  system.  The  same  results  do  not  exist  between  Lau(y)  and 
Lbu(y)  or  between  L  (y)  and  L.  (y)  if  an  optimum-spaced  array  is 
used.  av  Dv 
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